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EXPRESS MAIL NO. : EL387335666US 
ATTY. DOCKET NO. : SCHW-600 
PATENT 

SYSTEM AND METHOD FOR TAX SENSITIVE PORTFOLIO 

OPTIMIZATION 

COPYRIGHT NOTICE 

Material is contained in this application which is subject to 
copyright protection. The copyright owner has no objection to the 
facsimile reproduction of the patent disclosure by any person as it 
appears in the Patent and Trademark Office patent files or records, but 
otherwise reserves all rights to the copyrighted material. 

BACKGROUND OF THE INVENTION 

Recently, there has been a proliferation of software packages 
released that claim to help individuals make better investment decisions. 
Prior art investment software packages often utilize portfolio 
optimization methods to make specific security investment 
recommendations. Generally speaking, portfolio optimization refers to 
the process of maximizing the expected return on a portfolio for a pre- 
specified level of portfolio risk. Software packages that utilize portfolio 
optimization typically do not consider the tax consequences associated 
with investment decisions. As such, this class of portfolio optimization 
based software is only useful for making investment decisions within 
non-taxable accounts, such as 401(k) accounts. 

The importance of considering taxes in the portfolio optimization 
problem is well established in the prior art. Techniques for performing 
portfolio optimization when investors are subject to taxes are discussed 
in the prior art, but these techniques typically apply to investors making 
investment decisions within their taxable account only. Some 
techniques ignore existing investments for which capital gains taxes 
might be due if a sale is recommended. However, making investment 
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decisions within the context of the investor's entire portfolio, which can 
include both taxable and non-taxable investment accounts, is a generally 
accepted tenet of modern portfolio theory. 

What is needed is an investment system that makes investment 

5 recommendations within both taxable and non-taxable accounts and 
considers the consequences of those decisions within the context of the 
investor's entire portfolio. More precisely, an investor typically has a set 
of financial securities, an "investment universe/' that she can choose 
from when making investment decisions in her taxable account. This 

10 investment universe may or may not be the same as the universe of 
financial securities that she can choose from when making investment 
decisions in her non-taxable account. It is desirable to have an 
investment system that maximizes the investor's expected portfolio 
return after taxes for a level of portfolio risk that is determined by the 

15 investor's risk tolerance and provides investment recommendations 
within both taxable and non-taxable accounts. The output from this 
system could include specific dollar amounts, or investment 
percentages, to be invested in those financial securities included in the 
taxable investment universe and specific dollar amounts, or investment 

20 percentages, to be invested in those securities included in the non- 
taxable investment universe. 

SUMMARY OF THE INVENTION 
It is therefore an object of the present invention to provide an 
25 investment system that makes investment recommendations within both 
taxable and non-taxable accounts. 

It is another object of the present invention to consider the 
consequences of investment recommendations within both taxable and 
non-taxable accounts in the context of an investor's entire portfolio. 
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It is another object of the present invention to have an 
investment system that maximizes the investor's after-tax expected 
portfolio return for a level of portfolio risk that is determined by the 
investor's risk tolerance and provides investment recommendations 
within both taxable and non-taxable accounts. 

In accordance with one aspect of the present invention, a method 
is provided for characterizing an investment portfolio which includes the 
steps of: receiving data for one or more taxable investments; receiving 
data for one or more non-taxable investments; receiving investor profile 
information; performing an optimization which includes the data for the 
taxable investments, the data for the non-taxable investments and the 
investor profile information; and providing an investment 
recommendation. The investment recommendation may be, for 
example, with respect to asset classes or individual securities. The 
"investor profile information 7 ' could include, for example, an investor risk 
preference, a time horizon or tax bracket information. The optimization 
could be performed, for example, using a non-linear optimization 
routine. 

In accordance with another aspect of the present invention, a 
method is provided for characterizing an investment portfolio, 
comprising the steps of identifying a universe of taxable investments 
and non-taxable investments; receiving customer profile information; 
performing an optimization which includes the data for the taxable 
investments, the data for the non-taxable investments and the customer 
profile information; and displaying an investment recommendation. 

In accordance with yet another aspect of the present invention, a 
method is provided for characterizing an investment portfolio in a 
system in which a user communicates with a server by way of a 
computing device, comprising the steps of: identifying a universe of 
taxable investments and non-taxable investments; receiving customer 
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profile information; performing an optimization on the data for the 
taxable investments, the data for the non-taxable investments and the 
customer profile information in the server; and displaying an investment 
recommendation on the computing device. 

In accordance with still another aspect of the present invention, a 
method is provided for characterizing an investment portfolio in a 
system which includes a computing device, comprising the steps of: 
identifying a universe of taxable investments and non-taxable 
investments; receiving customer profile information; performing an 
optimization on the data for the taxable investments, the data for the 
non-taxable investments and the customer profile information in the 
computing device; and displaying an investment recommendation on the 
computing device. 

These and other objects, features and advantages will become 
apparent when considered with reference to the following description 
and the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 
FIGS. 1A through 1C are a series of exemplary questions for 

constructing an investor profile. 

FIGS. 2 through 2E are a series of asset allocation pie charts 

which correspond to the results of the investor profiles of FIGS. 1A 

through 1C. 

FIG. 3 is an overview of an apparatus for providing one 
embodiment of the present invention. 

FIG. 4 indicates a high-level application environment according to 
one embodiment of the present invention. 

FIG. 5 is a flow chart which illustrates the method according to 
one embodiment of the present invention. 
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FIG. 6 is a flow chart which indicates a first phase of the 
optimization process according to one embodiment of the present 
invention. 

FIG. 7 is a flow chart which indicates a second phase of the 
optimization process according to one embodiment of the present 
invention. 

FIGS. 8A and 8B are flow charts which indicate a third phase of 
the optimization process according to one embodiment of the present 
invention. 

FIG. 9 is a flow chart which indicates a fourth phase of the 
optimization process according to one embodiment of the present 
invention. 

FIG. 10 is a flow chart which depicts an overview of the data 
input, forecaster and optimizer features of the present invention. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
Referring now to Fig. 10, an embodiment of the present invention 
will now be described in greater detail. In accordance with the present 
invention, several functional areas are employed. Data are gathered in 
block 110, including investor data 112, data regarding non-taxable fund 
universe 114 and data regarding taxable fund universe 116. The 
investor data might include risk data, the customer's horizon data {e.g., 
years to retirement, etc.) and tax information {e.g., tax bracket, etc.). 

The Forecaster block 118 optionally performs an analysis of the 
expected performance of specified investments over a specified horizon 
and under specified conditions. 

The Optimizer block 120 uses the data from the input data block, 
specifies conditions to the Forecaster block 118 to obtain performance 
information about investment vehicles of interest, and then evaluates 
and optimizes an investment portfolio for the customer. The end result 
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of the method and system of the present invention is a recommendation 
of an investment portfolio that best satisfies the goals of the customer, 
taking into account the tax effects on the portfolio in view of the 
customer's tax status. 

The above functional blocks will be described in greater detail in 
the following sections in the context of investment vehicles of Charles 
Schwab and Co., Inc., ("Schwab") the assignee of the subject 
application. It is to be understood that the present invention is equally 
applicable to other investment vehicles and to those of other companies 
and entities, the following description being merely illustrative and not 
exhaustive of the features and details of the present invention. For 
example, where Schwab's mark Optimizer is referenced, one skilled in 
the art will readily understand that the following description is applicable 
to any such portfolio optimizer, however denominated. 

Inputs to the Optimizer 

An overview of the types of data that are input to the Optimizer is 
set forth in the following paragraphs. 

User-Defined Inputs 

Although the present invention need not involve Schwab 
accounts, it has been implemented at Schwab and therefore 
contemplates that a customer or client (who may or may not be the 
user) may have Schwab accounts and non-Schwab accounts. 
Information from all accounts should be input {e.g., number of shares 
and number of bonds). For each account the tax status is indicated 
{i.e., the user indicates whether the securities of the account in question 
are held in taxable or non-taxable accounts). 401(k) and 403(b) 
accounts should be identified. According to this embodiment, the 
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average cost basis for taxable positions should be indicated. It should 
be understood that average cost basis is just one way to capture cost 
basis information. Other accounting conventions such as first-in-first- 
out or specific identification could be easily accommodated. 

The client's initial cash position and the amount of cash to be 
added to each account during a specified time period, e.g., each year - 
taxable, non-taxable, and 401(k) - should be indicated. In addition, the 
client's federal and state tax rates should be indicated. 

Retirement profile information from the client should also be 
input, including date of birth, retirement age and retirement annual 
income goal from these accounts {i.e., excluding social security and 
other income sources). 

In this embodiment, a customer's risk tolerance, or risk level, 
must be input in some fashion. A customer's risk tolerance could be 
determined by filling out an investor risk profile such as the one set 
forth on FIGS. 1A through 1C. As noted on FIG. 1C, a customer's risk 
tolerance data may be combined with time horizon data to determine a 
corresponding asset allocation plan, such as the asset allocation plans 
set forth on FIGS. 2 through 2E. Alternatively, the customer could 
indicate her risk level by selecting one of a series of hypothetical 
investment plans, such as the asset allocation plans of FIGS. 2 through 
2E. 

Non-User Defined Inputs 

In addition to the current price of each security in the portfolio, 
there are seven input data files from which the Optimizer extracts data. 
The following paragraphs will discuss these data and their construction. 
Note that constructing these data involves a great deal of pre- 
processing. As we discuss below, this pre-processing involves executing 
programs on a periodic basis, e.g., quarterly or annually. The Optimizer 
simply extracts data from the data files that are the result of the pre- 
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processing, so none of the pre-processing needs to occur within the 
product itself. 

FIG. 3 illustrates one configuration of a system for providing data 
for the Optimizer. In this embodiment, customer account data, investor 

5 experience data and data regarding non-Schwab holdings are managed 
by Corporate Offering Data Services. The SEntry servers receive data 
from, and provide data to, both the IOS Web Server Complex and the 
Retail Web Server Complex. The Common Data Segment ("CDS") is a 
database (here, an Oracle database) in the web server environment 

10 which provides equity and mutual fund data to the IOS Web Server 
Complex and the Retail Web Server Complex. Branch representatives 
may submit requests to the Retail Web Server Complex, for example 
URL requests, and will receive data from the Retail Web Server 
Complex, for example via HTML. Alternatively, clients may submit 

15 requests to, and receive data from, the IOS Web Server Complex over 
the Internet. 

Mutual Fund Data File 
Return forecasts for individual securities are constructed from 
models of the way in which security returns move over time. These 
20 models depend on security type, i.e., mutual fund, stock, Treasury bond, 
cash. 

A future return forecast can be broken into three components. 
First, there is a measure of the expected abnormal performance - or 
expected alpha - on the asset. Second, there is a systematic 

25 component to return. This component depends on the returns to 
several common factors that drive the returns on stocks and bonds 
{e.g., the market risk premium, the market capitalization premium, the 
term structure premium, etc.) and the sensitivity of the security's return 
to movements in those common factors. Third, there is an unsystematic 

30 - or "surprise" - component to return. The unsystematic component is 
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that part of return that is completely unrelated to movements in the 
overall market {e.g., an unexpected labor strike at GM would likely have 
a negative effect on return that is unrelated to current market 
movements). Return surprises are unexpected, so estimates of this 
component of return are not needed for constructing expected return as 
is required to perform a portfolio optimization. 

In this section the details associated with estimating expected 
alpha, factor sensitivities and unsystematic returns for the universe of 
mutual funds that an investor might own are discussed. 
Estimating Expected alpha 
According to this embodiment, the annual expected alpha for 
mutual fund /is given by: 

annual E(a, ) =f 1 *.a i +f 2 * expense > , (1) 

where a, is the historical estimate of alpha for fund /; estimated 
over a prior period, e.g., the prior 36-month period, and expense, is the 
reported expense ratio. (The technique for estimating historical alpha is 
shown below.) The parameters f 1 and y 2 are the estimated return 
premiums for a one-unit exposure to a, and expense, , respectively. In 
this embodiment, the parameters f x and f 2 are estimated separately 
for domestic stock funds, international stock funds, and bond funds. 
(This estimation technique is also discussed below.) Expected alphas for 
international bond funds are computed using the domestic bond fund 
parameters. Parameters are not estimated separately for international 
bond funds because large samples of these funds are not available. 

According to this embodiment, the annual expected alpha for 
hybrid funds - funds that may own domestic stocks, international 
stocks, bonds or other securities - is determined using the following 
formula: 
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annual £fa,-) = (ff od(S x 3/ +ff ocks x expense ,)(% i/S stoofes,) + 
(f/* x a, + ff x expense ,){%nonUS stocks ,) + (2) 
(y! W x 5 ; +f| ow/ x expense ,)(1 - % i/5 sfoate, - %/7o/7 £/S stocks j) 

In equation (2), % £/S stocks, is the percentage of the fund's assets in 
U.S. common stocks and % nonUS stocks, is the percentage in foreign 
stocks. An implicit assumption in equation (2) is that the securities that 
5 are not classified as ordinary stocks by the data vendor, e.g., preferred 
stock and convertible securities, behave more like bonds than stocks. 

In this embodiment, equation (1) or (2) gives a fund's annual 
expected alpha for the first 2 forecast years. Following this initial 2-year 
period, 

10 

annual E(<Xj) = - expense, . (3) 

The reasons for doing this are twofold. First, it's very unlikely 
that a fund could maintain a positive alpha too long into the future due 

15 to competitive pressures. (In other words, a high positive alpha for an 
extended period is not consistent with a long-run equilibrium in the 
mutual fund industry.) In the absence of such an assumption, forecasts 
of future wealth will be extremely high for funds with high historical 
alphas. In addition, the Optimizer would heavily favor those funds with 

20 the highest alphas. Second, past research has shown that mutual fund 
returns are predictable for about a 1-year horizon. This favors giving 
high performing funds some credit for past performance, but not 
extrapolating this success too far into the future given the weak 
statistical relationship between past fund returns and fund returns more 

25 than 1 year in the future. 

The expected alpha computed from any of equations (1) - (3) is 

in terms of annual percentage return. The following formula is used to 

convert expected alpha to a monthly amount in decimal form: 
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monthly E (a,) = 

To distinguish among expected alpha estimates for the first 2 
years (24 months), define £ r i_ 2 4(a / ) as the monthly expected alpha for 
the first 24 months and E 25 _ H (a, ) as the monthly expected alpha for 
months 25 through the horizon month H. (In other words, H\s the 
number of months in the user's retirement horizon.) 

Estimation of f x and f 2 require a model of the way in which the 
returns on stocks and bonds behave over time. The modeling choice in 
the preferred embodiment is based on extensive academic research 
regarding the time series behavior of common stock and bond returns, 
especially Fama and French (1993). This research has shown that the 
following five-factor model does a good job explaining the returns on 
common stocks and bonds: 

r,,t ~ r f,t = a, + P/(W - r fit ) + s t SML t + vyMG t + ^ (5) 
t/TERM t + djDEFf + e t 

(Since mutual funds hold stocks and bonds, this model is also useful for 
mutual funds.) Here r i t is the fund return in period t, r f t is the risk- 
free return, a, is a historical measure of abnormal performance, 
r m t - r fit is the market factor, SML t is a small-minus-large-cap factor, 
VMG t is a value-minus-growth factor, TERM t is a factor that captures 
the return differential between long- and short-term bonds, and DEF t is 
a factor capturing the return differential between high- and low-grade 
bonds. Table 1 contains an exemplary set of factor data definitions. 
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Table 1: Exemplary Factor Data Definitions 



Variable 


Series 


r i 


US 30-davT-bill total return 


Domestic 
stock 




r m 


Large Cap Stock 


SML 


Large Cap Stock - Small Cap Stock 


VMG 


[(Large Lap vaiue — Large Lap urowuij -+- 

(biTiaii Lap vaiue omall v-,dp vjruwuiJJ/z 


International 
stock 




w 

' m 


World ex-US 


SML 


International Small Cap - World ex-US 


VMG 


International Value - International Growth 


Bonds 




TERM 


US Bond Horizon Premium 


DEE 


US Bond Default Premium 



At this point it is worth noting that different specifications can use 
different factors and/or different series to construct the factors. The 
5 important point is that the factors do a good job explaining the returns 
on the assets being modeled. In other embodiments, different factors 
may be required to model the returns on securities that are not modeled 
here. 

We use a two-step process to estimate f x and y 2 . For domestic and 
10 international stock funds, we estimate the following cross-sectional 
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ordinary least squares (OLS) regression for each quarter fin the sample 
period: 

n,t =lo,t+y i,f + y ^expense, + y v P/ + (6) 

where r /#f is the annualized return for the quarter, /Vis the number of 

funds in the sample for a particular quarter, the y kit 's k= 1,...,5 are the 

OLS regression slope coefficients (these are the coefficients being 
estimated) and e /#f is the regression error term. In this embodiment the 

regression (6) is estimated quarterly, but it could be estimated at other 
frequencies, e.g., monthly or annually. 

Prior to estimating (6), the coefficients s, and v, are 

estimated for each asset / using data over a prior period via the OLS 
regression (5), e.g., the prior 36 - 60 (at least 36 and no more than 60) 
month period. (The * symbol over a variable means that the variable 
has already been statistically estimated.) In this embodiment, the two 
bond factors are included in the equity fund time-series regressions 
because many stock funds own some bonds and because these factors 
do a good job explaining the return variation on these securities. By 
doing so, the coefficients on the stock market factors may be estimated 
with more precision. The average TERM and DEF premiums, however, 
are very small, so including the coefficients t, and d, in the cross- 
sectional regression (6) adds virtually no explanatory power to that 
regression. Including t, and d, in regression (6) would introduce 
estimation error - or statistical noise - into the estimation process, so 
these variables are not included. 

The cross-sectional regression (6) is estimated for every quarter 
in the sample period. The annualized time-series average, f k , of the 
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quarterly cross-sectional regression coefficients (the average of 
they\, 's) are then computed. 

In this embodiment, the same models are used to determine the 
expected alpha for international funds. International index data can be 
used in lieu of domestic index data when estimating regression (5). 
Unfortunately, international bond market data are not readily available. 
In one embodiment, domestic bond market factors are used instead of 
international factors when estimating sensitivities for international bond 
funds. 

Although the same approach can be used to estimate bond fund 
expected alphas, the models are different because different factors are 
better suited for pricing bonds than stocks. For each quarter in the 
sample, the following cross-sectional OLS regression model for bond 
funds is estimated: 



n.t=y o,f + Y i A + Y 2it expense, + y v p, 



According to this embodiment, the coefficients 3,,$/, t, and d, 
for each bond fund are estimated using data over the prior period via 
the following OLS regression: 

n, t - r fit = a, + p, {r mit - r fJt ) + t,TERM t + d,DEF t + e t . (8) 

The market risk premium is included in (8) because risky bonds 
tend to have some correlation with the overall stock market. The other 
stock market factors are not included in (8) because the estimated 
coefficients on these factors for bond funds tend to be extremely small 
and estimated with a great deal of statistical error. The cross-sectional 
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regression (7) is estimated using bond fund data for every quarter in the 
sample period, although other frequencies may be chosen. The 
annualized time-series average, f k , of the quarterly cross-sectional 
regression coefficients (the average of they kit 's) is then computed. 

Estimating Factor Sensitivities for Mutual Funds 

Since many mutual funds own both stocks and bonds, 
sensitivities to both the stock and bond market factors are estimated for 
the universe of mutual funds that a user may own. The methodology 
underlying the Forecaster and the Optimizer requires that index returns 
be available for a long time history. Unfortunately, historical return data 
on value and growth indexes prior to 1978 are not always readily 
available. For this reason, the VMG factor is not included in the return- 
generating models used by the Forecaster in this embodiment. Since 
return forecasts do not include a value-minus-growth factor, this factor 
is also excluded when estimating factor sensitivities. Since the market 
factor is correlated with the value-minus-growth factor, the regression 
coefficient on this factor - p - will capture some of the value-growth 
effect. Other embodiments could include a value-minus-growth factor 
or other factors that are considered important for explaining security 
returns. The following OLS regression model is estimated for each 
domestic mutual fund: 

r ir t ~r„ = */ +P/(W -r f , t ) + S,SML t + tfTERM + d,DEF + e t . (9) 

Because of data limitations (described below) the market 
capitalization factor SML is eliminated from the international mutual fund 
regressions in this embodiment. The following regression is estimated 
for each international fund: 
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r ht ~ r flt = a, + ^ orld -r fft ) + tfTERM + d,DEF + e t . (10) 

In this embodiment a world market index is used to measure the equity 
risk premium but U.S. bond market factors are used because bond 
market data are not readily available for international markets. 

The return data on individual securities for the mutual fund 
estimations are over a recent prior period, e.g., the most recent 18 to 60 
monthly returns (at least 18 months and no more than 60 months). For 
each fund, a,, p,- or pf° M , s, , i, and d, are estimated using (9) or 
(10). Note that for international funds s, =0 in this embodiment. Also, 
the estimate a, from the historical time-series regressions measures the 
historical abnormal performance on the fund and is one component of 
the expected alpha estimate, which is used to construct an estimate of 
future performance. 

A measure of unsystematic risk is required to simulate the 
unsystematic component of return. The standard deviation of the error 
term s e/l for each security /from regression (9) or (10) is the estimate 

of unsystematic risk. Most statistics packages report s e> , in their 

output. (Statistics packages refer to s eJ by different names. 

Sometimes s e i is called the standard error of the regression or root 

mean squared error.) 

To adjust for taxes on expected income and capital gains 
distributions, a tax multiplier is estimated- see equation (17). The 
inputs into the tax multiplier calculation are the pre-tax and after tax 
return on the fund over a prior period, e.g., 3 years. 

The mutual fund data file will be periodically updated, e.g., 

quarterly. The parameters of the expected alpha model should be 

periodically updated, e.g., once per year. 
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Table 2: Mutual Fund Data File 



Field 


Name 


i 


1 Inim 10 iHpnt*ifipr ftirWpr nr ni^in^ 

UMlv-JUC lUCllLIIICI ^LIL-lxGi \J\ \*\J^l\J J 


2 


Fund name 


3 


Fund category 


4 


^1-24 ( a /) 


5 


f25-//( a /) 


6 


p/ 


7 




8 




9 




10 


d, 


11 


S e,i 


12 


Tax multiplier 


13 


OneSource® designation (yes or no) 



5 Stock Data File 

The Optimizer will also input data on individual stocks. The stock 
data file is similar in format to the mutual fund data file, but the 
methodology for estimating model parameters is somewhat different. 
In this embodiment the expected alpha for individual stocks is 
10 assumed to be zero. In other words, past abnormal performance - 
performance in excess of a benchmark model - on stocks is not 
extrapolated into the future. This is done because there is very little 
evidence in the prior art indicating that abnormal performance on stocks 
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is predictable. However, it should be understood that an expected alpha 
model for stocks could be developed and incorporated here. 

The sensitivity of each stock's returns to the stock market factors 
are estimated via the following OLS regression over a recent period: 

r,.t ~ r f ,t = a, + P/ {r mrt -r f/t ) + s,SML t + e t . (11) 
International stocks are not included in the estimation in this 
embodiment, but given the data a similar procedure can be used to 
estimate parameters for these stocks. The standard deviation of the 
regression error term s ei from the estimation of equation (11) is the 

estimate of unsystematic risk for the stock. 

To adjust for taxes on expected income and capital gains 
distributions, a tax multiplier is estimated (see equation (17)). The tax 
multiplier has after-tax returns as one of its inputs. The inputs into the 
calculation of after tax returns are dividend yield and return over a prior 
period, e.g., the past 3 year period (see equation (20)). 

The stock fund data file should be periodically updated, e.g., 
quarterly. 



Table 3: Stock Data File 



Field 


Name 


1 


Unique identifier (ticker or cusip) 


2 


Stock Name 


3 


£/ 


4 


s, 


5 


S e,i 


6 


Dividend yield 


7 


Average annual historical return 
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Bond Data File 

The Optimizer will also input data on Treasury bonds. In this 
embodiment the expected alpha for individual bonds is assumed to be 
zero. So, like individual stocks, past abnormal performance on bonds is 
not extrapolated into the future. Again, an expected alpha model for 
Treasury bonds could be included here. 

The methodology chosen to model bond returns over time 
depends on the assumptions that are made regarding how the user 
rebalances their portfolio over time. For reasons discussed below, the 
preferred embodiment is monthly rebalancing, although other 
rebalancing frequencies can be accommodated. Furthermore, it is 
assumed that the user maintains the maturity of their bond portfolio 
over the retirement horizon at its current level. For example, a user 
who owns an intermediate-term bond is assumed to rebalance in such a 
way that she holds an intermediate-term bond in each period over her 
entire horizon {i.e., it is not assumed that the maturity decays over 
time.) 

As part of the pre-processing, factor sensitivities and the 
regression standard error are estimated for portfolios of short-, 
intermediate- and long-term T-bonds over a prior period, e.g., the 
previous 60 months. The following regression is estimated for each 
series: 

r i/t - r f/t = a, + 1, TERM t + d t DEF t +e t . (12) 

The standard deviation of the regression error term s eJ from the 

estimation of equation (12) is the estimate of unsystematic risk. 

To assign parameters to an individual T-bond, the T-bonds in the 
user's portfolio are classified as short-, intermediate- or long-term. In 
one embodiment, short-term bonds are defined as those bonds with less 
than 4 years to maturity. Bonds with maturities 4 to 10 years are 
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defined as intermediate term and bonds with 10 years to maturity and 
longer are classified as long-term. Obviously this is just one 
classification scheme and other schemes are valid. 

The bond fund data file should be periodically updated. 

5 

Table 4: Bond Data File 



Field 


Name 


1 


Asset Class (short-, intermediate, or 
long-term) 


2 




3 




4 


S e,, 



Cash 

Cash assets should also be included in the universe of securities 
10 from which the Optimizer makes recommendations. A logical choice for 
cash assets is money market funds. Money market funds have different 
returns due primarily to differences in fund expenses. The risk of 
money market funds and their covariances with other assets are roughly 
equal. Thus, in this embodiment a single money market fund is 
15 recommended. However, additional money market funds or other 
assets that proxy for cash investments can be included as additional 
securities. 

If the money market fund had no expenses, a reasonable 
estimate of return on this fund would be that it is equal to the return on 
20 the risk-free asset. Since it has an expense ratio, the pre-tax risk 
premium in this embodiment is estimated by the following: 



r sMMFj ~ r u = -e*pense SMMF , ( 13 ) 

-20- 

ATTORNEY DOCKET NO. SCHW-600 



PATENT 



where r SMMFit is the return on the Schwab Money Market fund. 

An estimate of the expense ratio on cash for customers with cash 
in their 401k account is required. Since data on expense ratios for 
5 money market funds are not readily available, equation (16) is used in 
this case as well. 

Factor Data File 

The methodology used to forecast future return scenarios 
10 requires a long history of monthly observations for the factors in the 
factor pricing models previously discussed. In addition to the model 
factors, the methodology developed later requires monthly observations 
on inflation and on a security that proxies for the risk-free asset. Table 
5 contains an exemplary factor data file. 



Table 5: Exemplary Factor Data File 



Variable 


Series 




US 30-day T-bill total return 


7T 


US inflation 


r m 


Large Cap Stock 


SML 


Large Cap Stock - Small Cap Stock 


VMG 


[(Large Cap Value - Large Cap Growth) + 
(Small Cap Value - Small Cap Growth)]/2 


w 

1 m 


World ex-US 


TERM 


US Bond Horizon Premium 


DEF 


US Bond Default Premium 
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Historical databases may have monthly observations on domestic 
series starting as early as 1926. Unfortunately, international stock 
return series are not typically available over such a long time frame. 
Instead of limiting the data to the more recent time periods where both 
domestic and international data are available, in this embodiment the 
U.S. bond market factors are used to model international bond returns, 
the SML factor is excluded for international funds, and the world market 
risk premium series is backfilled using the following approach: 

1. Using regression equation (14) below, estimate X 0 , 5^ and 
the standard deviation of the regression error term e, (s g ) using 
monthly data over the shorter time period where the world market risk 
premium series is available. 

-World r USTbill i r r S&P r US7W//\, c (14} 

2. For each month t beginning in 1/26, substitute the actual 
value for {r^y - rj^™*) and the estimates of £ 0 and X x from step 1 
into (14). 

3. Draw a random variable from a mean 0 standard deviation 
= 5 E normal distribution and substitute this variable for the term in 

equation (14). 

4. Repeat this process for each month during the period 
where world market risk premium data are not available. 

The main feature of the backfilling approach is that the backfilled 
international series is constructed so as to insure that the correlation 
structure between U.S. and international stocks observed in the actual 
data are preserved. The approach also links the mean and variation of 
the international equity risk premium to that of the U.S. While the 
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backfilling approach is no substitute for having real data, it's a good 
alternative to excluding ail data, including domestic data, for the earlier 
period where the international data are not available. More recent 
periods, e.g., the post-1970 period, are dominated by sub-periods of 
extreme performance and low volatility, relative to the long-run average. 



The factor covariance matrix Q (in this example, a 5x5 matrix) is 
constructed from the factor data file. It's used to construct the 
covariance matrix for the user's optimization problem. The T-bill and 
inflation series are excluded from the covariance matrix calculations. 
Using summation notation, the covariance between factor /and factory; 



where T is the total number of observations (T is the last 60 months), fj 



covariance matrix Q. 

The factor covariance matrix Q is more easily created using 
matrix algebra. Using matrix notation, 



Factor Covariance File 




T-l 



(15) 



is factor / = 1, 5, f ljt is the single element t= 1, Tin />, and /> 
is the sample mean of f, . When / = j, a f/j =aj , where af is the 
variance of factor /. Each o f /j is a single element in the 5x5 factor 



u = 



f'f -Tff 
T-l 



(16) 
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where f is the T x 5 factor matrix and f is the 5 x 1 vector of factor 
means. 

Recommended Fund Universe File 
There are two investment universes from which the Optimizer 
makes recommendations. The taxable universe includes those assets in 
the user's taxable accounts plus a Schwab determined mutual fund 
universe. Similarly, the non-taxable universe includes those assets in 
the user's non-taxable account plus a Schwab determined mutual fund 
universe. 

When choosing a fund universe from which to recommend funds, 
it's important to include funds from each asset class {e.g., large-cap, 
small-cap, international stock, bond, and cash) and funds that are 
expected to perform well in the future. Doing so allows the Optimizer to 
recommend broadly diversified portfolios of potentially high-performing 
funds. 

Prior research indicates that funds with good past performance 
and low expenses tended to perform well in the future. Schwab's 
mutual fund Select List™ consists of a small number of funds that have 
had either or both high past performance and low expenses in each of 
the major asset categories. Thus, the Select List™, or a subset of this 
list, is a good candidate to represent the fund universe. However, the 
Select List™ is just one example of a fund universe and any universe of 
funds could be used to practice this invention. 

It's important that the fund universe be relatively small. The 
optimization problem can be solved with a large number of funds, but to 
insure consistent and reliable advice it is preferable to have relatively 
small number of funds, e.g., 50 funds. The Select List™ typically 
contains 80 - 90 funds, so it is helpful to focus on a subset of the Select 
List™ universe. 
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As a starting point, all SchwabFunds® index funds 1 , OneSource® 
funds 2 and transaction fee funds with the exception of balanced funds, 
funds of funds, and any fund with a minimum investment amount above 
$3,000 are included. Balanced funds were eliminated because 
managers of these funds often make broad changes in their asset 
allocation. Funds of funds were eliminated because these funds do not 
have to report the underlying expense ratio on the funds they own. 
Only the expense layered on by the fund of funds itself is reported. 
Since the expected alpha model for mutual funds uses the expense ratio 
as an input, the expected alpha for these funds is measured with a great 
deal of error. Funds with minimum investment amounts above $3,000 
were eliminated because the optimizer will not be able to solve some 
cases because of investment minimums. By not allowing fund with high 
minimums into the universe, the number of cases where the optimizer 
fails because of a minimum investment requirement is minimized. 
SchwabFunds® that are printed in the publication but did not make the 

actual list are also excluded. 

This procedure may result in a universe size that is larger than 
preferred. The list can be paired down further in a variety of ways. One 
example would be to reduce the number of bond funds in the universe, 
because these funds tend to have similar risk/return characteristics. 

The fields required for the mutual fund data file are also required 
for the recommended fund universe file. In addition, the minimum 
investment amount for initial investments in a fund held in taxable and 
non-taxable accounts is required. See the description of the mutual 
fund data file for details regarding the estimation of model parameters. 



1 SchwabFunds® is Schwab's proprietary mutual fund family. 

2 Schwab's Mutual Fund OneSource® service includes over 1,000 no-load mutual funds 
from many of America's most prominent fund families all available without transaction 
fees through Schwab. 
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Table 6: Recommended Fund Universe File 



Field 


Name 


1 


Unique identifier (ticKer or cusip; 


2 


Fund Name 


3 


Fund category 


4 


^1-24 (« /) 


5 




6 


0/ 


7 


P World 


8 


5/ 
/ 


9 


f/ 


10 


^/ 


11 




12 


Tax Multiplier 


13 


Investment minimum taxable 


14 


Subsequent minimum investment 
taxable 


15 


Investment minimum non-taxable 


16 


Subsequent investment minimum non- 
taxable 


17 


OneSource® designation (yes or no) 



The cash fund, the Schwab Money Market fund, is also included 
5 in the recommended fund universe. In this case, the expected alphas 
are set to minus the expense ratio on the Schwab money market fund, 
the tax multiplier is not applicable (after-tax returns can be computed 
with knowledge of the user's tax rate) and the OneSource® designation 
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is yes. All other model parameters and the investment minimums are 
zero for cash. 

Model Portfolio Data File 
5 The Optimizer places constraints on the user's resultant asset 

allocation among stocks, bonds, and cash. In this embodiment the 
following data regarding Schwab's model portfolios are used: 



Table 7: Schwab's Model Portfolios 



Model 


%Large 


%Small 


%Int 


%Bonds 


%Cash 


Aggressive 


40 


25 


30 


0 


5 


Moderately 
Aggressive 


35 


20 


25 


15 


5 


Moderate 


30 


15 


15 


30 


10 


Moderately 
Conservative 


20 


10 


10 


45 


15 


Conservative 


15 


0 


5 


55 


25 



This is one example of desired asset allocations for different risk 
preferences. In other embodiments the number of models could be 
increased, the asset class list could be expanded, and/or the specific 
percentages could be changed. 

Tax Considerations 

The approach adopted in this embodiment to incorporating taxes 
into the portfolio problem is to consider the tax consequences of 
investment decisions. There are three tax sources that are considered 
in the investor's portfolio decision problem. First, investors may receive 
taxable dividend or interest income on their investments. Second, the 
tax consequences of the sale of securities held within a mutual fund are 
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passed on to fund shareholders. Third, when an investor sells a taxable 
security, capital gains taxes are incurred. To accurately consider these 
taxes one only needs to keep track of tax laws affecting the long- or 
short-term capital gains rate. The user will be responsible for specifying 
their ordinary income tax rate and whether the asset is held in a taxable 
or non-taxable account. 

Although it's possible to adjust variances and covariances for a 
tax effect, the preferred embodiment does not include such adjustments 
because they would have very little effect on variance-covariance 
estimates. Given the estimation error inherent in variance-covariance 
estimates, the added level of precision associated with adjusting these 
estimates for taxes seems like overkill. The existence of taxes would 
reduce variance calculations relative to their pre-tax levels. Moreover, 
the less tax efficient the asset the more the variance is reduced. 
Therefore, by not adjusting variances for tax effects a slight tendency 
towards the optimizer favoring less tax efficient assets is avoided. Other 
embodiments could include adjustments to the variance-covariance 
estimates to reflect tax effects. In the following, one embodiment for 
estimating after-tax returns is discussed. 
Mutual fund distributions 

There is a tax liability created from income or capital gains 
distributions that should be accounted for. These distributions take the 
form of dividends for stocks and interest for bonds. In addition to 
these, mutual funds also distribute short- and long-term capital gains. 

A popular measure of tax efficiency is the tax efficiency ratio: 

j. after -tax 

tax efficiency = ' x , 

where rf* e ~ ta * is the return on the asset prior to adjusting for taxes and 
r after-tax js the return a ft er adjusting for taxes. In a separate analysis, a 
tax efficiency measure was found to explain subsequent after-tax 
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returns to mutual fund shareholders. Given this, it seems logical that 
future after-tax return could be computed by multiplying future pre-tax 
return and the appropriate tax efficiency ratio. Unfortunately, this 
approach has its shortcomings. Most importantly, the tax efficiency 

5 ratio, as it's defined above, is only applicable when both pre-tax and 
after-tax returns are positive. If pre-tax returns are negative, for 
example, the tax efficiency ratio is greater than or equal to 1. A tax 
efficiency ratio greater than 1 implies that after-tax returns are bigger 
than pre-tax returns and this cannot be the case. Yet it is important to 

10 consider tax efficiency even when pre-tax returns are negative because 
the after-tax returns may be even more negative. 

In lieu of the tax efficiency ratio, a tax multiplier -is calculated as 
follows: 



15 77V, = —^r, < 17 ) 

(l + r^Y 12 



where the annual return lost to taxes for asset /; r[ axes , is 



a r pre-tax \ 
f . . _ ~ '/ >_ _ y 



20 



(l + zf*™) 



The denominator of the tax multiplier is taken to the 1/12 power 
to convert the tax multiplier to a monthly figure. The tax efficiency ratio 
and the tax multiplier are closely related. For assets where 
0 < tax efficiency < 1, the rank ordering of the assets in terms of tax 
25 efficiency is the same as the rank ordering using the tax multiplier in 
(17). The primary advantage of using tax multiplier instead of tax 
efficiency is that the tax multiplier is well defined even if the pre-tax 
and/or the after-tax return are negative. 
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Note that the tax multiplier is not a function of the user's tax rate. 
This is the case because data vendors cannot provide after-tax return 
estimates on a user-by-user basis. After-tax returns are typically 
computed by taxing income distributions at the maximum federal tax 

5 rate of 39.6% and capital gains distributions at 20%. These 
percentages have changed as tax laws have changed over time. 
Unfortunately, vendors who report after-tax returns do not always 
report the split between income and capital gains distributions. In this 
embodiment, a simple procedure to modify the tax multiplier to make it 

10 a function of the user's tax rate is illustrated. Other embodiments may 
exclude this step or use another technique for making the tax multiplier 
a function of the user's tax rate. 

Specifically, it is assumed that rf 3 * 65 = d,t where d, is the 
amount of the distribution for fund /and t is the tax rate. In this simple 

15 setting, d, represents the combined income and capital gain 

distributions and t is the weighted-average income and capital gains tax 
rate. To estimate d, given the after-tax return measure, an estimate for 

t is required. Given this, d, is given by: 



20 a,.V£±. (18) 



There are a variety of ways to estimate t However, the estimation of t 
should be done separately for stock and bond funds. 

Given a database with historical data on fund distributions, an 
25 average of the percentage of fund distributions attributable to capital 
gains and income could be computed for each year during the period in 
which after-tax returns are measured, e.g., each year during the prior 3- 
year period. The maximum federal income tax rate and the capital gains 
tax rate over the period during which after-tax returns are measured are 
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also recorded. In each year, t is estimated as the weighted average of 
the capital gains rate and the income tax rate where the weights are 
determined by the average of the percentage of fund distributions 
attributable to capital gains and income, respectively. The simple 

5 average of the yearly estimates i over the period in which after-tax 
returns are measured is the estimate of t used in (18). 

An alternative technique for estimating t is to estimate the 
historical split between dividend distributions and capital gains 
distributions using return and dividend yield information on a benchmark 

10 index. For example, the S&P 500 index could serve as the benchmark 
for stocks (or just large company stocks). An estimate of the 
percentage of total distributions attributable to income is the dividend 
yield divided by total return. One minus this amount is the estimate of 
the percentage attributable to capital gains. Given these estimates, the 

15 procedure outlined above can be used to obtain i . 
The adjusted tax multiplier is then 

TM?= ±— , (19) 

where t ol is the user's combined federal and state marginal tax rate on 
20 ordinary income. This is a conservative assumption in that by using the 
user's tax rate on ordinary income it is assumed that all distributions are 
either income oriented or short-term capital gains. An alternative would 
be to compute a weighted average income and capital gains rate for 
each user using the percentages of fund distributions attributable to 
25 capital gains and dividends determined from above. 
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r after-tax r pre-tax _ j. 
r i — '/ L oi x 



(20) 



where {— 



is the annual dividend yield on stock /and t oi is the user's 



tax rate on ordinary income. Given this, the tax multiplier for stocks is 
calculated via (17). 

After-tax return forecasts for stocks and mutual funds 
The expected return on each asset is input into the portfolio 
optimizer. To capture the likely impact of income and capital gains 
distributions over the life of the stock or fund, the following formula is 
used to compute the after-tax expected return for asset i: 

E after-tax ^ ] = ^ + £ pre-tax ^ ^ m _ _ ± _ ( 2 1) 

In the Forecaster, the pre-tax return on each security /also needs 
to be adjusted to account for taxes in month t: 



r after-tax 
r i,t 



= (l + r i f- t3X )xTM i -l. (22) 



When the asset in question is a mutual fund, the adjusted tax 
multiplier given by (19) is used in (21) and (22). 
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After-tax return forecasts for T-bond and cash returns 

In one embodiment, bond and cash returns are taxed at the 

user's ordinary income tax rate. In this case the expected after-tax 

return on bonds and cash is given by 

f jy. after -tax j = ^ pre-tax -j x (1 _ tgx ( 2 3) 

For cash, the tax rate is the user's combined federal and state ordinary 
income rate t oj . Treasury securities are exempt from state taxes, so 
the tax rate in this case is the user's federal tax rate t f . In the 
Forecaster, month t pre-tax returns on T-bonds and cash are adjusted 
for taxes as follows: 

a/ter-tax = pre-tax x (1 _ tgx rafe) ( 2 4) 

I ft I ft 

Selling a taxable security 

For asset /; the total taxes T, resulting from the sale of all or part 
of the user's current holdings are given by 

T, = V(w pi - w, )(1 - c, )t cg , (25) 

where V is the dollar value of the portfolio, w pi is the percentage of 
the portfolio invested in asset /prior to the sale, Wj is the "optimal" 
weight (or the weight after the sale takes place), <r,is the cost basis of 
asset /(expressed relative to the current price of the security) and t cg is 

the capital gains tax rate. 

These taxes reduce the current value of the user's portfolio. This 
reduction in portfolio value needs to be incorporated into the user's 
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optimization problem. A method to do so is shown in the section on 
building the Optimizer, below. 

BUILDING THE FORECASTER 

The goal of the Forecaster is to characterize the user's retirement 
age real wealth distribution (or alternatively the retirement income that 
this wealth could support). To do this, it is assumed that a handful of 
"factors" drive the performance of capital markets. In addition, security 
returns are modeled as a function of (1) the factor realizations {e.g., for 
the market factor, the actual return on a market portfolio of stocks 
minus the return on T-bills) (2) the asset's sensitivities to those factors, 

(3) a premium (positive or negative) for professional management 
(mutual funds only) and (4) a company- or fund-specific return surprise. 

The Forecaster extracts data from four input data sets: (1) the 
mutual fund data file, (2) the stock data file, (3) the bond data file and 

(4) the factor data file. To generate a single future return scenario for 

a security, historical factor realizations observed at a randomly selected 

date during the January 1926 to present time period, a random return 

surprise, the security's factor sensitivities and its expected alpha are 

substituted into the return-generating model. To preserve the 

correlation structure among stocks and funds, these same factor 

realizations are used to construct a corresponding future return scenario 

for each asset in the user's portfolio. A future portfolio return scenario 

and an associated real wealth level are then computed. (Future savings 

are also included in the real wealth calculation.) This process is 

repeated for each month in the user's retirement horizon, yielding a 

single observation for the distribution of retirement age real wealth. 

Repeating this entire process many {e.g., 5000) times creates a 

retirement-age real wealth distribution. This distribution is then 

characterized, e.g., by reporting its mean, median, and the 5 th and 95 th 
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percentiles. Also, statistics such as the probability that the user will 
meet their retirement objective can be computed. 

The following is a more detailed description of the methodology 
underlying the Forecaster, including justification for the modeling 
5 choices and an algorithm. 

Defining the return-generating model. 

In this section a method for using the time-series models (9) - 
(12) to generate future return scenarios is illustrated. Note that the 
models (9) - (12) are used to estimate the parameters p, , $f , s, , i, , 
10 d, and s e , . The Forecaster uses the following model to generate future 
pre-tax return forecasts for each asset in a portfolio: 

r iit =r ff +£(a / ) + p,(/v -^) + p7(/-^ -r^+^SMl, 
+ t,TBUi + djDEF. + e iit 

15 Note that the parameters p, , pf , s, , i if and d, are not all 

estimated for each asset. For example, pf is not estimated for 
domestic equity funds and p^ , t j f and d, are not estimated for 
individual stocks. Equation (26) can still be used to forecast returns as 
long as 0 is substituted for missing coefficients. For cash, all parameter 

20 estimates except £(a, )are set to zero. Note also that E(a, ) is 0 for all 

assets except mutual funds. 

For assets held within taxable accounts, the pre-tax return given 
by (26) is adjusted for taxes via (22) or (24). This calculation uses the 
tax multiplier. 

25 

Bootstrap the model 

The goal of the Forecaster is to characterize the user's retirement 

age real wealth distribution for a given portfolio. Of course the hope is 
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that the derived distribution adequately characterizes the "true" 
distribution of retirement-age expected real wealth. For this to be the 
case, the return generating model (26) must do a good job 
characterizing the expected returns and variances on stocks and bonds. 
In addition, the correlation coefficients among the estimated return 
series for the various assets need to be good estimates of the true 
correlations. 

The statistical properties - means, variances and correlations - of 
the factor return series used in (26) are the key drivers of those same 
properties for individual securities. The approach adopted in this 
embodiment for characterizing the retirement-age real wealth 
distribution hinges on the idea that the true distributional properties of 
the factor return series can be characterized by sampling from the long 
history of return data on the individual factors: statisticians refer to this 
technique as "bootstrapping." More precisely, historical factor 
realizations {i.e., actual factor returns) observed at a randomly selected 
date over the January 1926 to present time period are selected from the 
factor data file. Repeating this sampling technique many times insures 
that the statistical properties of the bootstrapped series are similar to 
their historical counterparts, e.g., the bootstrapped average equity risk 
premium is approximately equal to the average equity risk premium over 
the 1926 to present time period. 

To generate a future return scenario for a security, the 
bootstrapped factor return realizations from a single random draw, a 
security-specific random variate, and an estimate of expected alpha are 
substituted into (26). Expected alpha is set at f^Ca,) for the first 24 
months and f 2 s-//( a /) for the remaining months in the investor's 
retirement horizon. The security specific random return surprise is 
obtained by drawing a random variate from a mean 0 standard deviation 
= s er f normal distribution. 
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To preserve the correlation structure among the securities in the 
investor's portfolio, the same bootstrapped factor return realizations are 
used to compute a future return scenario for each asset in the portfolio. 
A future portfolio return scenario and an associated real wealth level are 
then calculated. This process is repeated for the number of months in 
the user's retirement horizon, yielding a single observation for the 
distribution of retirement age real wealth. Repeating the process many 
times generates a retirement-age real wealth distribution. 

There is evidence in finance literature that expected returns and 
variances vary over time. The methodology described thus far does not 
account for such variation. Properly accounting for time variation in the 
expected returns and variances of financial assets would have little 
effect on estimates of retirement-age return distribution statistics for 
users with long horizons. By not modeling this time variation, more 
precise estimates of retirement-age distribution statistics may actually 
be obtained because a signal-to-noise problem arising from the 
additional complexity brought about by modeling time variation is 
avoided. If shorter-term returns were being estimated or the path of 
the return series over time were considered important, then it might pay 
to accommodate time variation. Other embodiments might incorporate 
time-variation in expected returns and variances. 

Rebalancing 

The frequency with which to rebalance the portfolio must be 
determined. At one extreme one can choose not to rebalance, i.e., a 
buy-and-hold approach. Over long horizons, however, the portfolio 
return will be heavily influenced by assets with high expected returns 
{e.g., stocks), leading to extremely high estimates of future wealth. 
More frequent rebalancing mitigates this effect, so in this embodiment 
the portfolio is rebalanced monthly. 

-37- 

ATTORNEY DOCKET NO. SCHW-600 



PATENT 

When users rebalance they may incur transactions costs and 
taxes. Since these costs are difficult to estimate and tend to be small 
they are not included in future forecasts or in the optimization. 
(Transaction costs and taxes incurred on recommended trades are 
5 considered. Also, estimates of after-tax returns on mutual funds do 
include a tax on estimated capital gains distributions.) Having said this, 
other embodiments might use different rebalancing assumptions and 
incorporate, e.g., transaction fees on rebalancing trades. 

The effect of taxes associated with rebalancing on future real 
10 wealth forecasts and optimal portfolio weights can also be incorporated. 
However, the effects of taxes can be greatly reduced if the user: (1) 
uses future investment cash flows to rebalance; (2) offsets investments 
which increase in value with investments which decrease in value; 
and/or (3) does as much rebalancing as possible in non-taxable 
15 accounts. Thus, taxes associated with rebalancing are not modeled in 
this embodiment. 

The method for allocating future cash contributions to the assets 
in the portfolio must also be determined. If cash contributions were 
spread evenly across taxable, non-taxable and 401k accounts, annual 
20 contributions could be spread so as to maintain the original investment 
percentages of the portfolio. Unfortunately, annual contributions to 
401k accounts and other non-taxable accounts, e.g., IRA accounts, are 
capped. Thus, in many cases contributions cannot be spread so as to 
maintain the original investment percentages of the portfolio. Instead, 
25 contributions to an account are spread proportionally across the assets 
within each account to maintain the original allocation in that account. 
(Note that the forecasting logic described herein can also accommodate 
withdrawals - withdrawals are just negative contributions). 
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The Number of Simulated Paths 

The optimal number of simulated future return paths is a tradeoff 
between accurately estimating the true distribution of future wealth and 
run-time efficiency. With a small number of simulations {e.g., 100) the 
Forecaster will execute quickly, but the true distribution of future real 
wealth will be estimated with a great deal of error. This estimation error 
is most extreme in the tails of the distribution because there are only a 
small number of the simulated series that end up in the tails. For 
example, for a sample size of 100 there are only 5 observations in each 
5% tail of the distribution. Thus, it is preferable to simulate the model a 
large number of times. However, this is very expensive computationally. 

To help determine the number of simulated paths, experiments 
were performed using portfolios of differing security types. The 
objective was to determine the number of simulated paths required for 
the average ending wealth value of the simulated series to converge to 
its expected value and to have stable outcomes in the tails of the 
distribution. The number of simulated paths required depends on the 
type of assets held in the portfolio. A test group of stocks and mutual 
funds were selected. For each stock selected, 15 separate simulated 
distributions were estimated using a number of paths that ranged 
between 1,000 and 22,000. For each mutual fund, 15 separate 
simulated distributions were estimated using a number of paths that 
ranged between 1,000 and 15,000. A five-year horizon was used for 
each simulation. This horizon was chosen because the shorter the time 
horizon the greater the variability in terminal wealth, thereby giving a 
worst-case scenario. Most investors will have a much longer horizon, so 
the results for those horizon lengths will be even more stable and will 
converge more quickly. 

The analysis revealed that a stable solution for stocks was 
reached at 15,000 iterations; for domestic and international stock 
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mutual funds, a stable solution was reached at 5,000 iterations; for 
bond funds, a stable solution was reached at 1,000 iterations. 

In one embodiment 5,000 iterations is used. Another 
embodiment would allow the number of iterations to vary with the 
securities held in the portfolio. 

Converting Real Wealth to Annual Income 

The Forecaster estimates expected wealth, in today's dollars, at 
the user's horizon date. An estimate of the annual income that this 
portfolio could support during the user's retirement years is also 
reported. To avoid overwhelming the user with questions regarding 
their retirement investing and spending programs, a simplistic technique 
for converting real wealth into a lump sum is employed in this 
embodiment. Specifically, 4 percent of the real wealth level at the user's 
retirement age is the estimate of the retirement income this real wealth 
could support. The 4 percent number is based on a large body of 
practitioner literature indicating that a 4 - 6 percent withdrawal rate is 
appropriate for most investors. Other embodiments may vary the post- 
retirement withdrawal rate with the user's post-retirement horizon, 
portfolio risk, and the probability of running out of money before the 
user is expected to die. 

Calculating the Probability of Meeting the User's Goal 
For each of the 5000 future portfolio return scenarios, for 
example, there is an estimate of expected wealth and an associated 
income level in today's dollars. The distribution of annual income levels 
is sorted from lowest to highest. The annual income level nearest but 
not exceeding the investor's goal is then identified. The probability of 
meeting the user's goal is the rank position (l=low to 5000=high) of 
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this annual income level divided by 5000. If the stated goal is below 
rank 1, then the probability is 0. 

Asset Allocation 

The Forecaster will also display the asset allocation of the user's 
portfolio. The following describes one embodiment for the asset class 
definitions and mappings: 

Large Company Equity 

Mutual funds are placed in this asset class using two 
methodologies: 

1. Funds categorized as large- and mid-cap (includes blend, 
growth and value funds) are placed into this asset class. 

2. Specialty funds (communications, health, technology, etc.) 
with median market capitalization of $1 billion or greater are also placed 
into this asset class. 

Securities designated as domestic equity with market 
capitalization of $1 billion or greater are classified as large company 
equity. 

Small Company Equity 
Mutual funds are placed in this asset class using two 
methodologies: 

1. Funds categorized as small-cap (includes blend, growth 
and value funds) are placed into the asset class. 

2. Specialty funds (communications, health, technology, etc.) 
with median market capitalization less than $1 billion are also placed 
into this asset class. 

Securities designated as domestic equity with market 
capitalization less than $1 billion are classified as small company equity. 
International Equity 
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An exemplary list of categories for international equity funds are 
listed below: 



EXEMPLARY INTERNATIONAL EQUITY 
CATEGORIES 



DIVERSIFIED EMERGING MARKETS 
Diversified Pacific/Asia Stock 
Europe Stock 
Foreign Stock 
Japan Stock 
Latin America Stock 
Pacific/Asia ex-Japan Stock 
World Stock 



American Depository Receipts fADR's") are classified as 
international equity regardless of market capitalization. 
Fixed Income 

An exemplary list of categories for bond funds are listed below: 



EXEMPLARY FIXED INCOME 
CATEGORIES 

High Yield Bond 
Intermediate Government 
International Bond 
Long Government 
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LONG-TERM BOND 
Multisector Bond 
Muni National Intermediate 
Muni National Long 
Muni Short 
Muni Single State Intermediate 
Muni Single State Long 
Short Government 
Short-term Bond 
Emerging Markets Bond 
Municipal NY Intermediate 
Municipal New York Long Term 
Municipal California Intermediate 

Ultrashort Bond 
Municipal California Long Term 



Individual debt securities (e.g. corporate bonds, government 
bonds, etc.) are classified by debt type. See Table 8 for specific security 
categorization. 
5 Cash 

All funds held in cash and money market funds are placed in this 
asset class. 

Other 

Funds in this group might include 



EXEMPLARY OTHER CATEGORIES 

Convertibles 
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DOMESmC HYBRID 
International Hybrid 



In addition, individual items such as warrants, rights, options and 
convertible securities are also placed in this asset class. See Table 8 for 
specific security categorization. 



Table 8: Exemplary Asset-Class Categorization of Securities 



Security Type 


Asset Class 


Money Market Instruments 


Cash 


Commercial Paper 


Cash 


Bankers Acceptance 


Cash 


Certificate of Deposit 


Cash 


Money Market 


Cash 


Preferred Money Market/Muni Preferred 


Cash 


Money Market Fund 


Cash 


Reorganized Equity Security 


Equity 


When Issued Equity Security 


Equity 


User Defined Equity 


Equity 


Open Ended Mutual Fund 


Equity 


Equity 


Equity 


Real Estate Investment Trust 


Equity 


Stock 


Equity 


Common Stock 


Equity 


Preferred Stock 


Equity 


Convertible Preferred Stock 


Equity 


Cumulative Convertible Preferred Stock 


Equity 
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Corporate Bond 


Fixed Income 


Corporate Bond Coupon 


Fixed Income 


Corporate Bond Convertible 


Fixed Income 


Corporate Bond Original Issue Discount 


Fixed Income 


Mortgage Backed Security 


Fixed Income 


Government Agency Coupon 


Fixed Income 


Private Debt (No Secondary Market) 


Fixed Income 


Private Debt Convertible 


Fixed Income 


Government Agency 


Fixed Income 


Government Agency Original Issue Discount 


Fixed Income 


Government Agency Discount 


Fixed Income 


Collateralized Mortgage Obligation 


Fixed Income 


Mortgage Backed Pass-through 


Fixed Income 


Reorganized Debt Security 


Fixed Income 


When Issued Debt Security 


Fixed Income 


User Defined Debt 


Fixed Income 


Debt 


Fixed Income 


Asset Backed Security 


Fixed Income 


Municipal Bond 


Fixed Income 


Municipal Bond Coupon 


Fixed Income 


Municipal Bond Original Issue Discount 


Fixed Income 


Bond Fund 


Fixed Income 


Treasury (Government Debt) 


Fixed Income 


Treasury Bill 


Fixed Income 


Treasury Strip 


Fixed Income 


Treasury Coupon 


Fixed Income 


American Depository Receipt 


International 
Equity 


Option 


Other 
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Foreign Currency Option 


Other 


Index Option 


Other 


Broad Based Index Option 


Other 


Narrow Based Index Option 


Other 


Insurance 


Other 


Real Estate 


Other 


User Defined 


Other 


Index Master 


Other 


Unit Investment Trust 


Other 


OCC Pledged Security 


Other 


Basket 


Other 


Future 


Other 


Limited Partnership 


Other 


Derivative 


Other 


Commodity 


Other 


Right 


Other 


Warrant 


Other 


Closed End Mutual Fund 


Other 



Short positions are excluded from calculation of the asset 
allocation. 

Algorithm 

5 This algorithm is meant to serve as an example for how the 

Forecaster may be programmed. The algorithm implemented by one of 
skill in the art could vary, e.g., in order to increase the speed of the 
product. 

10 1. For each account (taxable, non-taxable and 401k), set V a t _ x = V a , 
where l^is the current total dollar value of the securities within an 
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account a. (Note that when the Forecaster is run for an optimal 
portfolio, V^ie must be reduced by the total dollar value of all capital 
gains taxes paid on recommended sell transactions and the total dollar 
value of all transactions costs.) 

2. For each account, set the weights invested in each asset to their 
initial value. When the forecast is being computed for the user's initial 
portfolio, the weights are set to those based on the user's original 
holdings. When the forecast is being computed for the optimal portfolio, 
the weights are set based on the optimal portfolio. 

3. For each month t = 1, retirement month, 

a. Randomly select a month £/from the 1/26 - present time 
period using a random number generator for a uniform distribution. 

b. For each account type: 

(1) For each security within an account type: 

(A) Substitute the factors for month u, including 
r fJU , into (26). 

(B) Substitute the security's estimates of the factor 
sensitivities into (26). If < 18 months of return data exist, then exclude 
that security. 

(C) For mutual funds, if t < 24 , then substitute 
(«/) int0 ( 26 )- Otherwise, substitute E 2S _ H (a,), into (26). For 

stocks and T-bonds, set £(a, ) = 0 . 

(D) Draw a random variable from a mean 0 
standard deviation = s e>i normal distribution. 

(E) Substitute this random variable for e l/t in (26). 
If the security is held in a taxable account, then adjust the pre-tax 
return using equation (22) or (24). 
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(2) Compute the portfolio return within an account: 
r .=7iv a ,r, „ , where r, , is the portfolio return on the account at 

time t, N= the number of securities in the account and w ai is the 
weight of security /in account a. 

(3) Adjust the return for inflation: 

real = 1 + r a,t wnere ^ js t h e month £/ inflation rate. 
a,t 1 + tc„ 

(4) Compute end-of-period real wealth: 

y f = V ait _iQ + cash contribution to that account. 

(5) Store V af . 



4. Add the end-of-horizon V a t 's. 

5. Repeat steps 1-4, e.g., 5,000 times. 

15 6. Compute distributional statistics for end-of-period wealth and 
retirement income: median, 5 th percentile and 95 th percentile. 
7. Compute the probability of meeting the goal as outlined above. 



THE OPTIMIZER 



In its simplest form, the portfolio optimization problem is to 
maximize expected return for a specific level of variance. Portfolio 
optimization is widely used in the pension fund industry where personal 
taxes are not considered. Different optimization routines are typically 
25 distinguished by their choice of modeling techniques for expected 

returns and the covariance matrix and by the constraints placed on the 
problem. The main distinguishing feature of the optimizer developed 
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herein is that an optimal portfolio is obtained for an investor who 
currently holds securities in both a taxable and non-taxable account. 

Methodology 

In this section the methodology for undertaking a tax sensitive 
portfolio optimization is displayed. The presence of taxes dramatically 
increases the complexity of the optimization problem, relative to the no- 
tax case. Specifically, the tax effect associated with capital gains or 
losses on securities is considered only in the case where securities are 
sold in a taxable account. Important implications of this are that the 
first derivative of the objective function is discontinuous; and that the 
objective function is a non-linear, and non-quadratic, function of the 
optimal weights. Moreover, the expected return on each asset depends 
on whether it's held in a taxable or non-taxable account. Because of the 
dual account nature of the problem, the optimizer must determine the 
allocation of funds within each account type, effectively doubling the 
number of parameters that must be estimated in the complete taxable 
case versus the case where no taxes are considered at all. 

Because the desire is for the solution to be tax sensitive and 
because a number of constraints are being imposed, nonlinear 
programming (NLP) will be required. The advantage of nonlinear 
optimization is that virtually any model can be specified, regardless of its 
functional form. There are two primary disadvantages: 

1. Nonlinear optimization is computationally expensive, since 
the functions involved may need to be evaluated hundreds or even 
thousands of times before the solution can be identified. The 
computational time required to solve the taxable case is also longer than 
that required to solve the no-tax case because there are twice as many 
parameters that need to be estimated in the taxable case. 
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2. Portfolio problems have very flat functions that are being 
optimized. In this situation the optimization will at times struggle to find 
an optimum. Thus, algorithms for what to do in the case where the 
optimizer stalls on a flat surface may need to be developed. It's also 
possible that the NLP routine will fail. Thus, this approach needs to be 
supported with error reporting mechanisms. 

In spite of these limitations, the optimization problem solved here 
is so general that this method should be adopted. To resolve the flat 
function problems identified in 2 above, an iterative solution algorithm is 
developed. This solution involves conditionally seeding the NLP routine 
with different starting conditions and, in the event that problems still 
exist, running the NLP again using information gathered from previous 
runs. An iterative procedure to incorporate investment minimums is also 
developed. The iterative procedure incorporates investment minimums 
without seriously hampering run-time efficiency. 

In the remainder of this section, a methodology for estimating 
expected pre-tax returns on the assets included in the optimization is 
developed first. Second, a technique for estimating the 
variance/covariance matrix without having to store the time-series of 
observations on each asset that might be included in the user's portfolio 
is shown. Third, the tax sensitive portfolio optimization problem is set 
up. Fourth, the iterative NLP solution technique is described. 
Estimating the expected pre-tax return on each individual asset 

An estimate of the pre-tax return for each asset in the portfolio is 
needed. The methodology for computing pre-tax expected return is 
illustrated using domestic stocks. 

The expected rate of return on each security in the user's 
portfolio is found by first applying the expectations operator to (9) - 
(12). In the case of domestic stocks, for example, applying the 
expectations operator to (11) yields 
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E(r, ) = E(r f ) + a,+ V,E{r m -r f ) + s,E(SML) . 

The unexpected return component e it in (11) disappears 
5 because the expected value of this variable is zero. Estimates of a, , 
p„ s n E(r f ),E{r mit -r fit ) and E{SML t ) are required to estimate 
expected return. The estimates of p,- and s, are the OLS estimates p,- 
and s, , respectively. The coefficient a, is meant to measure the 
expected abnormal performance on the security. This measure of 
10 abnormal performance is the security's E(a t ) , which is zero for stocks 
in this embodiment. The risk-free rate, market risk premium and the 
small-minus-large market capitalization premium are estimated using the 
long run time series average of these series 7> f r m - r f and SML, 
respectively. Making these substitutions yields the following formula for 
15 determining the pre-tax expected risk premium on individual domestic 
stocks: 

= ? f + V,r m -r f + SiSML, 

20 where r m - r f and SML are computed using the domestic versions of 
these series. 

For mutual funds an estimate of £(«, ) is required. As noted 
above, after year 2 E(a, ) is set to minus the expense ratio. [See the 
discussion of equation (3).] The optimizer needs a single input for 
25 F(a , ) . The following formula yields the average E(a , ) over the 
investor's retirement horizon H (in months): 

E{a,) = {[^^(a^f/^l + ^^a,)] 1 - 2 ^ Y" - 1 ' (27) 
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where f^Ca,) is the expected a over the first 24 months of the user's 
retirement horizon and E 25 _ H {aj) is the expected a over the remaining 
months. 

In general, the expected pre-tax return on any asset in this 
embodiment can be written as the following: 



E^ tax {r i ) = r f + -r, + $?r» -r f + s,SM+ (2g) 

ij TERM+ dj IJEF 



Note that for each security p, or will be zero. For many 
assets more than one factor sensitivity will be zero. Also, ^(a,) is non- 
zero for mutual funds only. 
Estimating the covariance matrix 

A method for estimating the covariance matrix without using the 
time-series of observations for each security is illustrated in this section. 
Given such a method, it is not necessary to store the time-series of 
observations for each of the assets that might enter the user's portfolio 
decision problem. 

Assume that a /C-factor model of the type defined in (9) - (12) 
describes the returns to all securities. [Model (9), for example, is a 4- 
factor model.] Let fibea Kx K factor covariance matrix, where 
<y f /j (the covariance between factor /and factor j) denotes the i,j th 

element of Q. When / = J, then o f/i =oj - the variance of the 
monthly returns on the I th factor. Furthermore, let s be the l\kN 
covariance matrix of the individual securities/positions included in the 
user's investment opportunity set. It is the covariance matrix s that is 
relevant for the user's optimization problem. 
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The individual elements u fJ in E can be constructed without the 
sample time-series using the /^-factor model. Define F to be the NxK 
matrix of factor sensitivities for N securities estimated from the K 
factors, fj , is the sub-vector of factor sensitivities for security /in F and 

f u a single element in fj {e.g., p, , §™ orld , s, , i, and d, ). The i,j th 
element of E is: 

where the final variance term equals zero except when i = j. Or, using 
matrix notation: 

Ofj ^=f i 'Of j +a|, (30) 

where the final variance term equals zero except when / = j. To 
construct the user-specific covariance matrix S in the example outlined 

in this embodiment, 

1. Define the factor covariance matrix Q to be a 5x5 (K=5) factor 
covariance matrix. Note that this covariance matrix includes both the 
U.S. and International stock market factor. For each security, either p ; 
= 0 or P/* 4 *** = 0 . For many assets more than one factor sensitivity will 
be zero. For example, an international mutual fund will have (3, = 0 
and Sj = 0 . 

2. For / = l,...yv, j = 1,...N, use (29) or (30) to compute the individual 
element in S. 

3. If /= j, the estimate sj #/ for security /is of . 
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Asset allocation constraints 

Asset allocation is a widely accepted tenet of modern finance and 
it is the foundation of Schwab's current investment advice philosophy. 
At the core of Schwab's asset allocation advice are the asset allocation 
percentages of the five model portfolios listed in Table 7. In the present 
embodiment the optimization problem is constrained to insure that each 
investor's resulting asset allocation is roughly consistent with one of the 
five model portfolios. The constraints do two things. Most importantly, 
they help to insure that recommended portfolios are diversified across 
asset categories. Thus, they help reduce the well-known tendency of 
optimizers to heavily recommend a single security with either an 
unusually high expected return or an unusually low variance. In 
addition, the constraints help to insure consistency with asset allocation 
models currently in existence. 

As an example, consider an investor who is classified as 
aggressive. The resultant asset allocation of the optimal portfolio is 
constrained to be roughly equal to the aggressive model portfolio 
allocations. More precisely, with the exception of cash the portfolio 
asset allocation is constrained to the percentages given in Table 7, plus 
or minus 5%. Cash is constrained to the corresponding asset allocation 
percentage plus or minus 1%. Since cash assets should be liquid, cash 
is forced into the taxable account to the extent possible. For example, if 
a customer has $10,000 of discretionary funds in her taxable account 
and the asset allocation constraint on cash implied that her overall cash 
position should be this amount or less, then her cash position would be 
held entirely within the taxable account. If the required cash position 
exceeded $10,000, then the remaining cash holdings would be held 
within the non-taxable account. 

In this embodiment, we do not make buy or sell 
recommendations of assets held within 401(k) or 403(b) accounts nor 
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do we make buy or sell recommendations for individual stocks or 
treasury bonds. These constraints, however, may be relaxed in other 
embodiments. 

To insure that there is a feasible space over which a solution can 
be obtained to each optimization problem, the asset allocation 
constraints apply to that component of the portfolio that the Optimizer 
can make recommendations. In other words, if the user holds 60% of 
her portfolio in individual stocks, for example, the asset allocation 
constraints apply to the remaining 40% of the portfolio. This is just one 
example of how to apply the asset allocation constraints in such a way 
as to insure a feasible space. 

The model portfolio concept is also helpful for communicating and 
specifying the user's risk tolerance. The user's optimization problem is 
solved for one of five different risk-return combinations based on the 
user's chosen attitude towards risk. These choices are labeled 
conservative, moderately conservative, moderate, moderately 
aggressive and aggressive following existing practice at Schwab. The 
decision to use five risk-return combinations and label them in this way 
is based on the following factors: 

• Risk is a difficult concept to convey, which makes the decision 
regarding labeling the user interface very important from a user 
experience perspective. Schwab's asset allocation models are the 
foundation of Schwab's existing advice platform. Customers and 
branch representatives are familiar with the model names: 
conservative, moderately conservative, moderate, moderately 
aggressive and aggressive. These names are also good names for 
characterizing the user's attitude towards risk. 

• The risk-return combinations on the efficient frontier - the set of 
portfolios offering the highest expected return for a given level of 
risk - are estimates. As such, they are subject to estimation error. 
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By having too many risk positions, there is a risk that the user will 
focus on differences in portfolio composition that are not meaningful. 
By keeping the number of risk positions small the user's attention is 
focused on portfolios that are more likely to truly reflect different 
expected risk-return combinations. 

The five risk-return combinations are determined by the asset 
allocation percentages corresponding to the five Schwab model 
portfolios. (See Table 7 for these percentages.) Therefore, the 
optimization problem can be solved for each model portfolio yielding five 
distinct sets of portfolio weights that represent alternate risk-return 
combinations. 

Setting up the optimization problem 

A key feature of a tax sensitive portfolio optimizer is that the 
expected return of each asset in the user's investment opportunity set 
depends on whether the asset is held in a taxable or non-taxable 
account. There are different approaches to setting up the optimization 
problem for this case. One embodiment is illustrated in this section. It 
will be apparent to those who have ordinary skill in the art that there are 
alternate ways to set up the optimization. What is common to these 
methods is that each requires an estimate of the amount of the asset to 
hold in the taxable account and the amount to hold in the non-taxable 
account. This implies that two parameters must be estimated for each 
asset included in the universe. The approach to this estimation problem 
adopted in this embodiment is to fold out the investor's universe. That 
is, each asset enters the universe once if it can be held in a taxable 
account and again if it can be held in a non-taxable account. 

In the remainder of this section, the specific steps taken to set up 
a tax sensitive portfolio optimization are outlined. 
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1. Combine multiple positions. Within each account - taxable, non- 
taxable and 401(k) - aggregate multiple positions in a single asset into 1 
position. For example, suppose there are two positions in mutual fund 
XYZ in a taxable account, 3 positions of that same fund in a non-taxable 
account, and 1 position in a 401k account. The two positions in the 
taxable account are rolled into 1 position by combining the shares in 
each position. Similarly, the 3 positions in the non-taxable account are 
combined into 1 position. Following this process, there are now 3 
positions in fund XYZ instead of 6 - one position within each account 
type. The primary advantage of combining assets into single positions is 
that run-time efficiency is improved because the number of parameters 
that must be estimated in the optimization process is kept to a 
minimum. 

2 . Determine the assets to be included in the user's optimization 
universe. 

• Within the taxable account, include: those assets that the 
customer holds in that account; those assets in the 
recommended fund universe file that customer does not 
already own; and the cash asset. 

• Within the non-taxable account, include: those assets that 
the customer holds in that account; those assets in 
recommended fund universe file, with the exception of 
municipal bond funds, that customer does not already own; 
and the cash asset. 

• Include those assets held within a 401k in the non-taxable 
universe. 

• Define /Vas the number of assets in the universe. 

Note that this method of folding out the universe implies that the 
number of parameters that the optimizer must estimate is atieasttwo 
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times the number of funds in the recommended fund universe. (An 
equivalent method for estimating investment percentages in taxable and 
non-taxable accounts is to include a second parameter measuring the 
percentage of the security to be held in the taxable/non-taxable account 
5 for each security that can be held in both types of accounts. The 

number of parameters to be estimated is the same under both of these 
methods.) 

3 . Group assets within each account by constraint type. 
10 • Group 1: The non-cash assets held in the taxable account 

where a buy or sell can be recommended. Define N t as the 
number of assets in this group. 

• Group 2: The non-cash assets held in the taxable where only 
a se//can be recommended. Define N x as the number of 

15 assets in this group. 

• Group 3: The cash asset(s) held in the taxable account. 
Define N 3 as the number of assets in this group. Note that 

N 2 = 1 in this embodiment. 

• Group 4: The assets held in the taxable account where 

20 neither a buy nor a sell can be recommended {e.g., stocks 

and T-bonds). Define /V 4 as the number of assets in this 

group. 

• Group 5: The non-cash assets held in the non-taxable 
account where a buy or sell can be recommended. Define 

25 N s as the number of assets in this group. 

• Group 6: The non-cash assets held in the non-taxable 
account where only a sell can be recommended. Define /V 6 

as the number of assets in this group. 
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• Group 7: The cash asset(s) held in the non-taxable account. 
Define N 1 as the number of assets in this group. Note that 
N 1 = 1 in this embodiment. 

• Group 8: The assets held in the non-taxable account where 
neither a buy nor a sell can be recommended {e.g., stocks, 
T-bonds, and 401k assets). Define N 8 as the number of 

assets in this group. 

8 

Note that W = £/V/ . 

/=i 

4. Construct the recommended asset universe. The 8 groups listed 
above combine to form the recommended universe of assets. For 
purposes of exposition, the universe is arranged in such a way that the 
group 1 positions are at the top and group 8 positions are on the 
bottom. Obviously, this sorting method is only one way in which a 
mapping between the recommend asset universe and optimal 
investment amounts can be created. 

5 . Dynamically determine cost basis equation for each asset held in 
a taxable account. To do this, sort lots in security /held in the taxable 
account from highest cost basis to lowest. 



Lot 


Cost 
basis 


Current # of 
shares 


Current 
Price 


1 


q (highest) 




Pi 


























L 


c L (lowest) 


s L 
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x, = 



V[\r\ax(w pi -w,,0)] 



where V\s the value of the portfolio, w pi is the weight of asset /in the 
original portfolio, and w, is the optimal weight. The cost basis can then 
be written as a function of the optimal weights: 



Cj = minCSi ,x / )-c 1 + min[s 2 , max(;r , - s t ,0)] • c 2 + ■ ■ ■ 
+ m\n[s L , max(X/ -s 1 -s 2 s L ,0)] • c L 



(31) 



In this embodiment the position with the highest cost basis should be 
sold first. Other embodiments may consider alternate conventions such 
as first-in-first-out. 

6. Determine the asset allocation boundaries. Define the 6x1 vector 
b of asset allocation percentages for the model portfolio corresponding 
to the user's chosen level of risk: 



A' 



where: th corresponds to the percent allocated to asset class 1, defined 
as large cap; bz corresponds to the percent allocated to asset class 2, 
defined as small-cap; th corresponds to the percent allocated to asset 
class 3, defined as international; £ 4 corresponds to the percent allocated 
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to asset class 4, defined as bonds; £ 4 corresponds to the percent 
allocated to asset class 5, defined as cash; and th corresponds to the 
percent allocated to asset class 6, defined as other. Note that b 6 = 0 
for each model portfolio. These percentages are taken from Table 7. 
Also, define the vector of plus-minus percentages as 



Pi 



Pe 



In this embodiment p x = 1% and all other values for pare set to 
5%. Note that the values in p imply that the asset allocation 
percentages are allowed to vary by plus or minus 5% for all asset 
classes except cash, which is only allowed to deviate by plus or minus 
1%. 

7 . Define the weight vector and determine initial investment 
percentages. The vector of optimal investment percentages - or 
weights - corresponding to the user's asset universe is represented by 
the following: 



w 



w 



N 



, , , % amount invested in position i . 
(Note that w, = • ) 
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The vector of initial investment percentages is similarly defined as 



w p = , 

where w pj is the user's current investment percentage in position /. 

8 . Determine discretionary investment percentages in taxable and 
non-taxable accounts. For group 4 and group 8 assets the optimization 
is constrained so that no buy or sell recommendation is made. This 
means that in many cases the optimizer will have discretion over less 
than 100% of the total portfolio. If N x + N 2 + l\/ 3 > 0 , then the percent 
of discretionary money in the taxable account is defined as 

Dtax = Yu w P i '< otherwise, D tax = 0. Similarly, if N s + N e + N 7 > 0 , 

i=i 

the percent of discretionary money in the non taxable is defined as 

/+/V 5 +/V 6 +/V 7 -1 

D non-tax = X" > otherwise = 0. 

/=/Vi+/V 2 +/V 3 +/V 4 +l 

9 . Determine the upper and lower boundaries on the portfolio 
weights. In traditional optimization problems each weight on the N 
assets in the universe is typically allowed to vary between zero and one. 
This insures that only long positions in the asset are considered. The 
upper and lower boundaries on the portfolio weights can be further 
restricted in this problem, improving the run-time efficiency of the 
product. The weights within the 8 groups of assets are constrained as 
follows: 
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Group 1: For each w, , i = 1, ■ -,N X , 0 < w, < D tax . The 
weights on group 1 assets are constrained so they can not 
exceed the discretionary amount in the taxable account. 

2 

Group 2: For each w,, i = \ + N 1 , 2 + N t , — , ^ N < > 

i=\ 

0 < w, < w pl . This set of constraints guarantees that the 
optimizer will not recommend a buy on a group 2 asset. 

2 2 3 

Group 3: For each w,, i = 1 + ^N, ,2 + ]T /V, ■ -^N, , 

/=i 1=1 /=i 

0 < w, < C UB , where the cash upper boundary C UB is set as 
follows: C UB = {b x + p^Pte + D non _ tax ). (One could set 
the cash upper boundary to 

Cug = min{(^ + A )(D tax + D non _ tax ), D tax } if it is more 
efficient to do so.) The cash upper boundary is set to take 
advantage of the fact that the model portfolio percentages 
can be used to constrain the weight even further than would 
be the case without these constraints. The upper boundary 
calculation includes non-taxable cash because cash is being 
forced into the taxable account whenever possible. 

3 3 4 

Group 4: For each w lf i = l + £/V,,2 + 2]/V / , ■-,Y J N i > 

/=i /=i /=i 

w, = w pi . This constraint insures that a group 4 position will 
remain fixed in the user's portfolio. 

4 4 5 

Group 5: For each w,, / = 1 + J^N,- ,2 + , • -,2^/ ' 

/=i /=i /=i 

0 < w, < D non _ tax . The weights on group 5 assets are 
constrained so they can not exceed the discretionary amount 
in the non-taxable account. 
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5 5 6 

• Group 6: For each w,, i = l + X yv /' 2 + £ /v /' "'IX ' 

/=i /=i /=i 

0 < < w p/ . This set of constraints guarantees that the 

optimizer will not recommend a buy on a group 6 asset. 

• Group 7: Cash is placed in the taxable account whenever it is 
possible to do so. To accomplish this, the upper boundary 
on non-taxable cash is determined as a function of the upper 
boundary on taxable cash. Define the upper boundary on 
non-taxable cash to be C UB inon -tax ■ For each w t , 

/ = l + £ n, , 2 + X N, , ■ ■ ■, X N, , , 0 < w, < C m non^ , where 
i=i /=i /=i 

Cub, non-tax casn =0 If C UB * O m , i.e., the upper boundary 
on total cash is less than or equal to the percent of 
discretionary funds in the taxable account. Otherwise, 
C UB , non-tax = C UB - D Tax , i.e., the upper boundary on non- 
taxable cash is equal to the upper boundary on total cash 
minus the amount of discretionary funds in the taxable 
account. 

7 7 8 

• Group 8: For each w,, i = l + ^W / ,2 + £/V / , ■■■ l ^ N i > 

i=l /=1 /=1 

w, = w pi . This constraint insures that group 8 positions will 
remain fixed in the user's portfolio. 

It is possible to pare down the upper boundaries on some of the 
assets even more by taking advantage of what is known about the asset 
allocation constraints. 
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10. Create the matrix of asset allocation weights. Define the N x 6 
matrix a containing the asset class weights of each of the N assets in 
each asset class / as 



a = 



7 /VI 



'16 



'A/6 



If /is a group 4 or group 8 position, the individual elements 
a ffl j = I,- -,6 in a are set to 0. This implies that the style of each of 
these assets is not considered in the asset allocation constraint. If /is 
10 not in group 4 or group 8, the individual element a iJf j = 1,- -,6 in a is 

set to 1 if the asset belongs to asset class J and 0 for all other elements, 
The 6x1 vector a'w contains the asset class weights for each of the 6 

asset classes. 

15 11. Calculate pre-tax expected return. For each asset N i n the fund 
universe, determine E pre ' tax {r,-^- the pre-tax expected rate of return - 
using (28). 

12. Calculate returns after taxes on distributions. For group 1-4 
20 positions, calculate E pre ~ cap - ga/n [r,] - the return after adjusting for 

expected taxes on subsequent distributions but before capital gains 
taxes - using (21) or (23). 

13. Set up expected return vector. Define the /Vxl vector of 
25 expected returns as 
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ppre-cap. gain -j 



r - 



E ore-cap. gam r r i 
l//V 1 +/l/ 2 +/V 3 +/V 4 J 
c pre-tax r r n 
^ l//V 1 +/V 2 +/V 3 +/V 4 +l J 



p pre-tax r r l 
c L'/Vi +/V 2 +/V 3 +/V 4 +1 J 



14. Sef matrix to determine capital gains/losses. When an asset 
that the user currently owns is sold a capital gain or loss is incurred. To 
capture this tax effect first define 



max{w p2 -w 2r 0) 

A w = 

msAw pN - w Nl Qi) 

\, is an Nx\ vector where each element takes on the value w pi - w> 

when a sale of that security takes place in the taxable account or 0 
otherwise. In addition, define the N*l vector 1-c that contains 1 minus 

the cost basis for each position: 



l-c = 



1-q 
1-c, 



1-c 



N 
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When asset /is not held or it's held in a non-taxable account, set c, = 1 . 

The total dollar tax Tincurred due to selling securities in the 
user's portfolio is given by 

T = t cg x min{l/[(l - c)'A w ], 0} , (32) 

where t cg is the capital gains tax rate. In this embodiment capital losses 
on assets are allowed only to the extent that these losses in aggregate 
are offset by capital gains. That is, the amount of income that can be 
offset is set to zero. The tax laws currently allow for a net loss of 
$3,000 to be offset against ordinary income in a given year. Moreover, 
additional losses can be carried forward. A user's ability to take 
advantage of these losses for tax purposes is influenced, among other 
things, by the level of taxable income. Since it is difficult to estimate the 
benefits of selling assets at a net loss, the conservative assumption that 
no benefits are received is adopted in this embodiment. This assumption 
has the effect of discouraging the sale of assets for the purpose of 
realizing a capital loss. Another embodiment would be to cap the net 
amount at $3,000 instead of 0. The drawback of this is that the user 
may use the product within the same year and (32) would again allow a 
$3,000 deduction. 

15. Account for transaction costs. At Schwab, transaction fees are 
levied on mutual fund buys/sells for funds not in the mutual fund 
OneSource® program. Expected portfolio return should be adjusted for 
these fees, commissions on stocks, and commissions on Treasury bonds. 
Transaction fees for non-OneSource® funds at Schwab are as follows: 
■ 0.7% of assets for amounts below $15,000, with a $39 
minimum. 
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■ 0.2% of assets for amounts at or above $15,000 plus 0.7% 
for the first $14,999, not to exceed a max of $149. 
Unfortunately, precisely modeling this pricing structure within the 
optimization is extremely cumbersome because a fund can be held in 
multiple accounts and fund purchases within each account are a 
separate transaction. This means that the size of the transaction fee 
depends on the purchase amount in each account. Modeling this would 
have a deleterious affect on product performance. To avoid this, in this 
embodiment the transaction fee for each mutual fund / is modeled as 

gf nd = max[min(rf, x V Q x \w pi - 1 x 0.007, d, x 149), d, x 39] , (33) 

where stakes on the value of 1 if the fund is non a OneSource® fund 
and 0 if it is a OneSource® fund. Note that this formula is applicable 
for fund sales as well as fund purchases. 

Stock commissions at Schwab vary based on the medium used to 
make the trade - phone, branch or web - and the number of trades an 
investor makes in a month. Since the preferred embodiment of this 
product is a web-based tool and is designed for long-term buy-and-hold 
type investing, the standard web commission fees of $29.95 per trade or 
$0.03 per share for trades over 1000 shares is used. Define the total 
number of shares for a stock position /as S, . Commissions on stocks 
are modeled as follows: 

gf** = max(S, x \w pi -w,\x 0.03,29.95) , (34) 

Commissions on Treasury bonds at Schwab are a flat $49 per 
transaction. Therefore, if position /is a T-bond, then 
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gT-bond _ 4 g # (35) 

Given these definitions of g, , the total dollar value of all 
transaction fees are 

5 

G = tfl, 06) 

where g is an Nx 1 vector containing the g, 's and 1 is an Nx 1 vector 
of ones. 

10 

16. Calculate after-tax portfolio return. First, recognize that the 
expected portfolio return pre-capital gain taxes e 9ain [r p ] is 

^pre-cap gain^^ _ w ^ _ 

15 

The expected portfolio return after all taxes E^-^Vp^ 



E after ~ tax [r p ] = 



i 

V-T-G 



H 

-1 



V-T-G\h 



V 

i 

(l + w'r)-l, 

J 
20 

where H\s the total number of months in the user's investment horizon. 

17. Compute the variance-covariance matrix. Follow the steps 
outlined in the previous section to compute the Ak/Vcovariance matrix 
25 S. 
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According to one embodiment, the optimization is formulated as a 
single-period problem. Those skilled in the art will recognize that 
another embodiment would be to formulate the problem as a multiple- 
period problem. The latter approach may be desirable in the case 
5 where risk and expected returns on securities are allowed to vary over 
time. The objective of the optimization is to determine the weight w, 
to be invested in each of the N positions so that the following function is 
minimized: 

10 min -?.w'Sw 

w 2 

subject to 

(b - PX^te* + & non-tax) < a'W < (b + p)(D tax + 

^non-tax ) 

15 w'l = 1 

/Vi+yVj+^+yv, /v 1 +/v 2 +/v 3 +/v 4 

/=i /=i 

and the Group 1 through Group 8 constraints identified in step 9, 

where X - the risk aversion parameter - is set by the mathematical 
20 programmer. As set forth in more detail below, x, is determined by the 

customer's risk level. 

The term w'Zw in (37) represents the matrix calculations 
required to compute portfolio variance. When formulating the user's 
optimization problem, it is standard to multiply the variance term by the 
25 constant Vz. This simplifies the formulas for the first and second 

derivatives of the objective function without having any affect on the 



'V-T-G 
V 



(l + w'r)-l 



(37) 



-70- 

ATTORNEY DOCKET NO. SCHW-600 



PATENT 

optimal solution. If the inputs into (37) are in terms of percentages, 
however, this constant should be 0.005. 

The term in brackets is the portfolio expected return, adjusted to 
account for taxes and transaction costs. The objective function is 

5 formulated to capture the idea that the user dislikes variance and likes 
expected return. 

The first constraint is the asset allocation constraint, applied to 
only that portion of the portfolio that the Optimizer has the discretion to 
recommend trades. The second constraint makes certain that the 

10 weights sum to one. The third constraint restricts the percentage 

amount held in taxable accounts to its original value. This insures that 
the amount of money in the user's non-taxable account does not change 
and thereby rules out the possibility of triggering other tax ramifications. 
The remaining constraints are the group 1 through group 8 constraints 

15 identified in step 9. In summary, the constraints of this embodiment 
rule out short selling, fix the weight on those assets (such as individual 
stocks and assets in a 401k account) for which buys or sells are not 
recommended, and limits the investment in those assets not included in 
the Schwab determined fund universe file to no more than their original 

20 weight. 

Although the risk of the user's portfolio is primarily determined by 
the asset allocation constraint, a value for the risk aversion parameter x 
must be specified to solve (37). The following model is used to choose k 
for each optimization run: 

25 

\n(X) = 3.33 + 0.65 x t oj - 2.3 x \u{risk level) , (38) 

where the risk level in this embodiment varies from 1 (conservative) to 5 
(aggressive). These risk levels correspond to the asset allocations 
30 depicted in FIGS. 2A through 2E, as described above. Note that 
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equation (38) gives the natural logarithm of X. In order to determine X, 

one calculates the exponential of this value. 

The model for X was determined by experimenting with a 
prototype. Specifically, for an aggressive investor with rollover cash and 
a 35% tax rate the value for X that roughly yielded a portfolio with the 
highest Sharpe ratio was determined. For lower levels of risk 
(moderately aggressive, moderate, etc.), the value for X that yielded an 
asset allocation that roughly corresponds to the asset allocation 
percentages for the model portfolio associated with that risk level was 
determined. 

For each risk level, the tax rate was lowered by 10% and the 
value for X that yielded a portfolio variance similar to the portfolio 
variance at the corresponding risk level for the 35% tax rate was found. 
This process continued for successively lower tax rates. An OLS 
regression was fit to the resultant data yielding the results presented in 
(38). (The r-square value for this regression is 0.997.) 

The optimization problem (37) should be solved using a non- 
linear optimization routine. An NLP solution technique is discussed next. 

An Iterative NLP Solution 

Several factors heavily influenced the design of the solution 
algorithm. The most important of these factors are as follows: 

1. The objective function is non-linear and non-quadratic. 

2. The constraint set includes inequality constraints. 

3. All constraints are linear. 

4. The first derivative of the objective function is discontinuous. 

5. The variance-covariance matrix for the full sample is singular or 
near-singular (because the number of assets is near the number of 
observations used to estimate the input parameters and because there 
are redundant assets in most cases). 
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6. Mutual funds have minimum investment requirements. 

These factors imply that the traditional quadratic programming 
approach to solving the portfolio optimization problem will not work. 
Instead, an iterative NLP solution is used in this embodiment. Since the 

5 variance-covariance matrix is nearly singular, it is often difficult to 

identify the maximum - or minimum - of the objective function when all 
assets are evaluated. The iterative solution algorithm is designed to 
solve these difficult cases. An iterative solution algorithm is also used to 
accommodate minimum investment amounts. 

10 In this embodiment the NLP algorithm is a quasi-Newton 

algorithm where the first and second derivatives are approximated. A 
method that considers the second derivative is used because the 
feasible solution set may be non-convex. (A convex set is a set where a 
line drawn between any two points in the feasible space is contained 

15 within the set.) While evaluating the second derivative is important 
when the feasible solution space is non-convex, no NLP algorithm can 
guarantee that a global maximum/minimum has been found. This last 
point implies that one should be careful when communicating the 
optimization process to the user. 

20 Note that traditional portfolio optimization problems that do not 

consider capital gains taxes or investment minimums are quadratic, 
implying that the feasible space is convex. This class of problems can 
be solved using methods that look at the first derivative only. By 
evaluating the second derivative, the computations necessary to find a 

25 solution are increased considerably. The matrix of second derivatives - 

or Hessian - has ~ ^ distinct elements, each of which has to be 

N 

approximated at each iteration of the NLP (the NLP may require as few 
as 10 iterations or as many as thousands of iterations to find a solution). 
This large increase in the amount of computations will negatively 
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influence the speed of the routine relative to an efficiently coded routine 
for the simpler case. 

The fact that analytical approximations to the derivatives are 
used means that a solution can be obtained despite the discontinuity in 
the first derivative. At points very near the point of discontinuity, the 
analytical approximation may be a poor one. An alternative 
embodiment is to use a non-derivative NLP algorithm (e.g., grid search) 
instead of a quasi-Newton algorithm. Non-derivative methods tend to 
be slow and it is often difficult to identify that a maximum/minimum was 
found. For this reason and because quasi-Newton methods are well 
established in the optimization literature and in practice, the preferred 
embodiment of the present invention uses a quasi-Newton approach. 
There are a variety of quasi-Newton routines that one could implement, 
which are well known in the art. 

The iterative solution algorithm can be separated into four distinct 
phases. In phase 1 an attempt is made to resolve flat function problems 
by starting the NLP routine with different sets of initial values. The 
preferred embodiment uses three different sets of initial values. This 
can be reduced to a single set to improve run-time efficiency. If the 
function does not optimize with these values, in phase 2 the solution for 
the best case from phase 1 is used as the initial values and the 
optimization is rerun including only those assets with nonzero weights. 
Given an optimal solution from phase 1 or phase 2, in phase 3 the 
optimization is rerun to account for the minimum investment values. In 
phase 4, the solution is reported to the user. 

In the flowcharts referenced in the following discussion, the 
following notation is used to represent solution status: STAT=1 indicates 
an optimal solution has been identified; STAT=2 indicates a feasible 
solution can not be identified; and STAT=65 indicates the attempted 
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stepsize is less than the precision {i.e., the flatness of the function 
precludes a definite solution). 

The first phase is illustrated on FIG. 6. Initial values (weights) 
are first set to zero. If STAT=1, the optimization proceeds to phase 3. 
If not, the initial values are set to 1. Once again, If STAT=1, the 
optimization proceeds to phase 3. If not, the initial values are set to 
1/N. Again, If STAT=1, the optimization proceeds to phase 3. If not, 
the optimization proceeds to phase 2. 

Phase 2 is depicted on FIG. 7. In phase 2, it is first determined 
whether any of the solutions from phase 1 are feasible. Then, it is 
determined which of the various solutions from phase 1 produces the 
highest value for the objective function. Then, the NLP is run on the 
subset of assets with nonzero weights from the best phase 1 solution. 
The optimization then proceeds to phase 3. 

Phase 3 is depicted on FIGS. 8A and 8B. Phase 3 uses a series of 
branches to identify a solution that incorporates the minimum 
investment required for each asset (to fully account for this would 
require a mixed-integer NLP, which would likely require far too much 
computational time to solve). This phase receives as input an optimal 
(STAT=1) solution from phase 1 or 2. In the first case, 3.0, all weights 
in the optimal solution exceed the minimum so the optimization routine 
proceeds to phase 4 for cleanup and reporting. 

If at least one weight in each asset class meets the minimum 
investment requirement, then the optimization proceeds to phases 3.1. 
and 3.2. In phase 3.1, the NLP is run on the existing solution, setting 
the lower bound for all nonzero weights equal to their minimum 
investment value. In phase 3.2, the NLP is run on the existing solution 
where all assets that have weights less than the minimum investment 
value have their upper bound clamped to zero. The solutions for phases 
3.1. and 3.2. are compared and the optimal solution with the highest 

-75- 

ATTORNEY DOCKET NO. SCHW-600 



PATENT 

value for the objective function is selected. This provides the final 
solution for phase 4. 

As noted in FIG. 8B, if at least one asset class has no weights 
meeting the minimum investment, the NLP is rerun in phase 3.3 as 
follows. For all asset classes with at least one asset meeting the 
minimum investment, those assets meeting the minimum requirement 
are included, their lower bound is set to the minimum and those 
securities not meeting the minimum are excluded. For asset classes 
with no assets meeting the minimum, the one with the highest weight is 
identified and its minimum value is set to the investment minimum. All 
other assets in this class are excluded. If the NLP solves correctly, the 
optimization proceeds to phase 4; otherwise, the optimization proceeds 
to phase 3.4. 

In phase 3.4, positive and negative deviations from the asset 
class bounds are included in the optimization. The sum of the 
deviations is included as a penalty in the objective function whose 
impact is varied by the parameter Gamma, which is set to 100 in this 
embodiment. Negative deviations in the cash targets are penalized at 
2x in this embodiment. Variables are initialized as in case 3.3. This 
method generally assures a solution when the minimum investment 
constraint creates problems compounded by the inability to move cash 
and first stage solutions from taxable to nontaxable accounts. In order 
to preclude "lumpy" solutions, assets with high minimum investments 
(above $3,000) are excluded from the recommended fund universe in 
this embodiment. The maximum allowable deviations from the initial 
targets are currently unconstrained. 

Phase 4 is illustrated on FIG. 9. At this point an optimal solution 
has been identified from the NLP. Phase 4 is a cleanup and reporting 
phase which makes the transformations necessary to link the results of 
the optimization process with the final user interface. There are three 
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identifiable cases to be addressed where a solution can not be reported 
to the user: 

1) The NLP can not solve due to the lack of funds— it is possible, a 
priori, to determine when there is not a feasible solution due to the 
fund minimums and how the discretionary funds are allocated across 
the tax/nontax accounts. A reporting mechanism should be designed 
to account for this case. 

2) The NLP does not indicate a solution in any phase before the 
minimum investment criteria are considered— given the robustness of 
the iterative solution method, this is not likely to occur but it cannot 
be ruled out. Thus, an error reporting mechanism must provide for 
the case where the NLP just does not work. 

3) There are situations in phase 3 where it is possible for the NLP to 
fail. Although the flexibility of the 3.4 technique should preclude 
most all of these, an error is still possible. An error reporting 
mechanism is needed for this possibility. In these cases, however, 
modifying the solution approach may resolve the error. 

Special Cases 

As a result of the criteria used to eliminate funds from the solution 
set when attempting to satisfy the minimum investment constraint, it is 
possible to define an infeasible problem. There are two cases currently 
addressed: 

1) If the user has initial cash in the nontaxable account and the 

preliminary solution has amounts invested in the fund universe that 
do not meet the minimum requirements, the funds in phase 3 will be 
clamped to zero and the cash constraint will preclude holding these 
funds in nontaxable cash. Thus, in this case, there is nowhere to 
invest the nontaxable discretionary funds. The current embodiment 
addresses this problem by looking at the upper bound for fund 
universe assets in the folded-out space. If the sum of the upper 



-77- 

ATTORNEY DOCKET NO. SCHW-600 



PATENT 

bounds is less than the proportion of discretionary nontaxable funds, 
the upper bound on nontaxable cash is set equal to the original 
proportion of cash in the nontaxable account. 
2) Similarly, in both the taxable and nontaxable account, if new assets 
are obtained by selling all of an existing asset, then the weight of the 
asset sold is clamped to zero. If the recommended weights were 
such that they did not meet the minimum investment requirement, 
then phase 3 of the program, in some cases, will clamp their 
maximum to zero. Once again, this potentially leaves no place to 
invest the funds within a particular account. To solve this, the sum 
of the upper bounds for existing mutual funds and the funds in the 
Schwab determined fund universe is compared with the discretionary 
funds available in that account. If the sum of the upper bounds is 
less then the discretionary funds for the account, then the upper 
bound on all previously owned assets is reset to its original value 
(allowing funds to revert to their original position if they are 
insufficient to be moved to a new asset). 

Implementation and User Interface 

The Optimizer could be implemented in a variety of ways. For 
example, Schwab representatives could input customer data, perform 
the optimization at the PC or server level, and then provide the 
information to the customer. Alternatively, the Optimizer and all 
necessary non-customer data could be provided to the customer's PC, 
and the customer could perform the optimization. In the preferred 
embodiment, the Optimizer application software resides on a server, for 
example an IOS Web Server as shown in FIGS. 3 and 4. 

A representative "clickflow" diagram, indicating how a user could 
navigate one implementation of the Optimizer, is shown on FIG. 5. 
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Starting from a planning overview screen, the user is welcomed. First- 
time users are prompted to enter portfolio information, to select 
accounts, to indicate an annual contribution level and retirement 
information and to enter risk data. The user is then guided to a profile 
5 summary screen. Users who have already entered the foregoing 

information would proceed directly to the profile summary screen. The 
user would then (optionally) proceed to one or more portfolio value 
forecasting screens. 

The user could then choose to optimize her portfolio. The user 
10 could view asset allocation comparisons, change assumptions subject to 
certain constraints, and would then have the option of re-optimizing her 
portfolio. The user's account could be charged for using the Optimizer, 
either on a per-session or a per-optimization basis. 

One skilled in the art will appreciate that various modifications or 
15 applied examples are conceivable which are within the scope of this 
invention. Accordingly, the scope of this invention is not limited to the 
previously-described embodiments. 
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WE CLAIM: 

1. A method for characterizing an investment portfolio, comprising 
the steps of: 

receiving data for taxable investments; 
5 receiving data for non-taxable investments; 

receiving investor profile information; 
performing an optimization which includes the data for the 
taxable investments, the data for the non-taxable investments and the 
investor profile information and which takes into account capital gains or 
10 losses on taxable investments which would be sold; and 
providing an investment recommendation. 

2. The method of claim 1, wherein the step of providing an 
investment recommendation includes the step of providing weights 
which indicate a level of investment for asset classes of the taxable 

15 investments and non-taxable investments. 

3. The method of claim 1, wherein the step of providing an 
investment recommendation includes the step of providing weights 
which indicate a level of investment for securities of the taxable 
investments and non-taxable investments. 

20 4. The method of claim 1, wherein the step of providing an 
investment recommendation includes the step of providing dollar 
amounts to be invested in asset classes of the taxable and non-taxable 
investments. 

5.. The method of claim 1, wherein the step of providing an 
25 investment recommendation includes the step of providing dollar 
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amounts to be invested in securities of the taxable and non-taxable 
investments. 

6. The method of claim 1, wherein the investor profile information 
comprises an investor risk preference. 

7. The method of claim 1, wherein the investor profile information 
comprises tax bracket information. 

8. The method of claim 1, wherein the investor profile information 
comprises time horizon information. 

9. The method of claim 1, wherein the investment recommendation 
comprises a first recommendation for asset classes of the taxable 
investments and a second recommendation for asset classes of the non- 
taxable investments. 

10. The method of claim 1, wherein the investment recommendation 
comprises a first recommendation for securities of the taxable 
investments and a second recommendation for securities of the non- 
taxable investments. 

11. The method of claim 1, wherein the step of performing an 
optimization includes the step of performing numerical differentiation. 

12. The method of claim 1, wherein the step of performing an 
optimization includes the step of employing a quadratic optimization 
routine to obtain a solution in a non-quadratic environment. 

13. The method of claim 1, wherein the step of performing an 
optimization comprises an iterative non-linear optimization routine. 
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14. The method of claim 10, wherein the first recommendation for a 
taxable investment has at least one asset class in common with the 
second recommendation for a non-taxable investment. 

5 

15. The method of claim 10, wherein the first recommendation for a 
taxable investment has at least one security in common with the second 
recommendation for a non-taxable investment. 

16. The method of claim 13, wherein the optimization routine 

10 comprises a quasi-Newton algorithm wherein a first derivative and a 
second derivative are approximated. 

17. The method of claim 13, wherein the optimization routine 
comprises a first subroutine of attempting to resolve a flat function 
problem by running the routine with different sets of initial values. 

15 

18. The method of claim 17, wherein a function does not optimize 
with any of the sets of initial values used in the initial step, further 
comprising a second subroutine of: 

taking a solution for a best case; and 
20 re-running the optimization routine including only those 

investments with nonzero weights. 

19. The method of claim 17, wherein an optimal solution is found in 
the first subroutine, further comprising a third subroutine of re-running 

25 the optimization routine to account for minimum investment values. 

20. The method of claim 18, wherein an optimal solution is found in 
the second subroutine, further comprising a third subroutine of re- 
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running the optimization routine to account for minimum investment 
values. 

21. A method for characterizing an investment portfolio, comprising 
the steps of 

5 identifying a universe of taxable investments and non-taxable 

investments; 

receiving customer profile information; 
performing an optimization on the data for the taxable 
investments and for the non-taxable investments which takes into 
10 account capital gains or losses on taxable investments which would be 
sold; and 

displaying an investment recommendation. 

22. A method of characterizing an investment portfolio in a system in 
which a user communicates with a server by way of a computing device, 

15 comprising the steps of: 

identifying a universe of taxable investments and non-taxable 

investments; 

receiving customer profile information; 

performing an optimization on the data for the taxable 
20 investments and for the non-taxable investments in the server which 
takes into account capital gains or losses on taxable investments which 
would be sold; and 

displaying an investment recommendation on the computing 

device. 

25 23. A method of characterizing an investment portfolio, comprising 
the steps of: 

communicating with a server by way of a computing device, 
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identifying a universe of taxable investments and non-taxable 
investments; 

receiving customer profile information; 

performing an optimization in the server on the data for the 
taxable investments and for the non-taxable investments which takes 
into account capital gains or losses on taxable investments which would 
be sold; and 

displaying an investment recommendation on the computing 

device. 

24. A method of characterizing an investment portfolio with a 
computing device, comprising the steps of: 

identifying a universe of taxable investments and non-taxable 
investments; 

receiving customer profile information; 

performing an optimization on the data for the taxable 
investments and for the non-taxable investments in the computing 
device which takes into account capital gains or losses on taxable 
investments which would be sold; and 

displaying an investment recommendation on the computing 

device. 

25. The method of claim 22, wherein the receiving step comprises 
receiving the customer profile information at the computing device. 

26. The method of claim 24, wherein the receiving step comprises 
receiving the customer profile information at the computing device. 
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ABSTRACT OF THE DISCLOSURE 

A method and apparatus are provided which account for tax 
impact in predicting investment return. Data for taxable investments, 
data for non-taxable investments and investor data are received and an 
optimization is performed which includes all of these data. An 
investment recommendation is provided to the investor based on the 
optimization. 
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Determine what type of investor you are 
and match a plan to your goals. 

The type of plan that best suits you depends on the type of investor you are — your 
Personal Investor Profile. Are you an aggressive investor, a conservative investor or 
somewhere in between? We Ve designed this questionnaire to help you find out quickly 
and easily. By answering these questions, you can combine 1) the time period you want 
to consider for investing, your "time horizon;' and 2) how you feel about the risks of 
investing, yourJYisk tolerance',' to define your investor profile. 

Once you decide which Personal Investor Profile best describes you, it's easy to choose the 
best way to allocate your assets and diversify your portfolio to meet your goals. 
The questionnaire assumes that the money you plan to invest will be used for a long-term 
purpose such as retirement or a single purpose such as saving for a home. Read it over 
and check the answers that apply to you. Then follow the instructions to select an asset 
allocation plan. 

I plan to begin withdrawing money from my investments for major 
needs within: 

Less than 3 years 1 point 

3-5 years 3 points 

6-10 years 7 points 

i 1 years or more ? 10 points 

——points 

^) Once I begin withdrawing funds from my investments, I plan to spend 
all of the funds within: 

Less than 2 years 0 points 

2-5 years 1 point 

6-10 years 4 points 

1 1 years or more 8 points 

points 

Subtotal A 

Add your points from questions 1 and 2 and enter here __ 

If your score is less than 3, STOP HERE. 

A score of less than 3 indicates a very short investment time horizon. For such a short 
time horizon, a relatively low-risk portfolio of 40% short-term (average maturity of five 
years or less) bonds or bond funds and 60% cash is recommended, as sfock investments 
may be significantly more volatile in the short term. Go directly to the short-term asset 
allocation chart. 



If your score is greater than 3, please continue. 

FIG. 1A 



€> 



I would describe my knowledge of investments as: 

None 0 points 

Limited . . 2 points 

Good . 4 pomts 

Extensive . . .... . . , 6 points 



. points 



When I decide how to invest my money, I am: 

Most concerned about the possibility of my 

investment losing value . 0 pomts 

Equally concerned about the possibility of my 

investment losing or gaining value 4 points 

Most concerned about the possibility of my 

investment gaming value . 8 pomts 



, points 



© 



Review the foliowing list and select the investments you currently own or have 
owned in the past. Then choose the one with the highest number of points and 
enter that number. 

(Example You Money market funds or cash equivalents 0 points 

, now ^ st ° ck Bonds and/or bond funds 3 points 

funds and have in r 

the past purchased Stocks and/or stock funds 6 points 

international secu- International securities and/or international funds . 8 points 

nttes. Your score x M , , 

would be 8) (Maximum possible score = 8) points 

Consider this scenario: 

Imagine that over the past three months, the overall stock market lost 25% of its value. 
An individual stock investment you own also lost 25% of its value. What would you do? 

I would 

Sell ail of my shares 0 points 

Sell some of my shares 2 points 

Do nothing 5 points 

Buy more shares 8 points 

___ points 

^} Review the following chart: 

We've outlined the most likely best and worst case annual returns of five hypothetical 
investment plans. Which range of possible outcomes is most acceptable to you or best suits 
your investment philosophy? 

(The figures are hypothetical and do not represent the performance of any particular investment.) 
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Investment plans 


Avg. annualized 
return (1 year) 


Best case 
scenario { I year) 


Worst case 
scenario (1 year) 


Points 


Investment Plan A 


12% 


16.3% 


-5.6% 


0 points 


Investment Plan B 


9 09 


25.07c 


-12.1% 


3 points 


Investment Plan C 


10 4% 


33 69c 


-18.2% 


6 points 


Investment Plan D 


11 19c 


42 8% 


-24.0% 


S points 


Investment Plan E 


12.59 


50.0% 


-28.2% 


10 points 



points 

Subtotal B 



Add your points for questions 3 through 7 and enter here 



Now you're ready to 
determine your Personal 
Investor Profile and 
choose a corresponding 
asset allocation plan. 

o 

Enter Subtotal A here 

This number represents your 
time horizon score 

Enter SuBlOFAL B here : 



This number represents your 

RISK TOLERANCE SCOre. 



Time Horizon Score 



3- 4 points I Spoints } 7-9 points 1 10-12 points! 14-18 points 

12 At G conservative "- *' ■ ; v,V:^M 

13 - \ .. : ^ " f" " ' ' v . . ^| 

is . 9 • - 

16 ^%\\}-p^7?V\ . - ■ : '/-"M 

17 ■ ; B-A/Idderately Canservatftft*^ 

18 '[-l[ : 9\ 9 . ■ , ' ;i 

19 l^W^'^J 9 '9 



GO 



C~ Moderate 



D- Moderately Aggressive 



Plot your time horizon score and 
your risk tolerance score on the 
Personal Investor Profile chart and 
locate their intersection point. Then 
find the corresponding asset allocation 
pian on the following page. 



E-Aggressive 



FIG. 1C 



Now find your corresponding asset allocation plan. 

The following allocation plans match the Personal Investor Profiles on the previous page— A through E. 
Now that you know your own profile, go to the corresponding chart below to determine how to allocate 
your assets Then go to Step 2 to choose a solution that suits your plan 

Short term Bonds AQ% 



For investors who want current 
income and a high degree of 
stability tn their investments 
Gfven 3 short-term time horizon, 
it is prudent to invest in asset 
classes that experience relatively 
minor price fluctuations, such as 
short- term fixed -income funds 
and money market funds 
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Cash 60% 



o 



Conservative 



Average annuaf return 
(1370-1996) 9 44% 
Best year 21 81% 
Worst year -1 25% 

For investors who want 
current income and stability 
and aren't concerned about 
increasing the value of 
their investments 




Cash 25% 



FIGi. 2A 



Bonds 55% 



Moderately 
Conservative 

Average annual return 

(1970-1996}: 1065% 

Best year 25 57% 

Worst year: -6.54% Lar 8 e " 
company 

For investors who want stocks 
current income and stability, 2C% 
with some increase m the 
value of their investments 



Bonds 45% 




FIGi. 2B 



Cash 15% 



FIGi. 2C 
FIGi. 2D 
FIGi. 2E 
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Moderate 



Average annua! return 
(1970-1996) 1155% 
Best year. 29.76% 
Worst yean -12 95% 

For long-term investors who 
don't need current income 
and want reasonable but 
stable growth Some fluctua- 
tions are tolerable, but they 
want less nsk than the overa! 1 
stock market 




Cash I0 ; ; 



Large-company 
stocks 30% 



©Moderately 
A . international 

Aggressive stocks 



Average annual return 
0970-1996) 1141% 
Bestyear. 36 63% 
Worst year -19 14% 

For long-term investors who 
want good growth and don't 
need current income A fair 
amount of risk is acceptable, 
but not as much as if they 
invested exclusively in stocks 



Bonds 15% 




Cash 



Urge-company 
stocks 35% 



|^ Aggressive 

Average annual return 
(1970-19961 1189% 
Best year 42-39% 
Worst year -23 82% 

For long-term investors who 
want high growth and don't 
need current income. 
Substantial year-to-year 
volatility in value k acceptable 
in exchange for a potentially 
high long-term return 



Cash 5% 




Small-company 
stocks 25% 



Large-company 
stocks 40% 



Clients 




Non-Schwab 
Holdings 



Investor 
Experience 
Planning 
Data 
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1.13 High-Level Application Environment 



Portfolio Optimizer 

High-Level Application Environment 




i —J JSEntry 

Libraries 

I Firewall #2 



y 

JSEntry Servers 



V 




DB2 
Database 
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FIG. 5 



3.3 Page "Click flow" Diagram 

The following diagram shows the click ftou for all pages used in Optimizer. 



Portfolio Optimizer 
- "Clickflow" 



Planning 
Overview 



Welcome 



I Welcome Back 



pes 

j Portfolio Selection 




VL 



Select Your 
Accounts 




OG5 Enter Ann! 
Contribution 



Retirement Info 
1 




;ocb your 
Contribution 



CB 




Risk Profile 


i - 


i 






■ Profile Summary 



Cs Risk Assess 
Questionnaire i 



T 



° Portfolio Value 
Forecast Calcs 



KEY: 

->- link 

► nav button 

C^nT) button 



xxx i "P arent " 
1 page 



xxx Ji pop-up or 
^ "link" page 



ttnked from 
Steps Bar 

linked from 
sub-Steps Bar 



T 



°Your Forecast - 
Current Portfolio 



TJ 



° Optimized 
Portfolio 




Asset Allocation 
Comparison / 



Optimize Change 
Assumptions ^ ^ 

w QoptimizeJ) 



Optimization 
Constraints , 



Plan Summary 



Phase 



1. Try \anous initial values If speed becomes an 
issue, the number of initial values attempted can be 
reduced and or the precision of the \LP can be 
reduced. 



1 



Set initial 
weights = 0 



N ST AT o* 



i .2 Set initial 



j \ST*T= 



\$T\T 



3 Set initial j nstat*. 
weights = i N. j 



[ 3 



NSTAT o\ 
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2. Are any of the solutions 
from Phase 1 feasible? 



False 



REPORT 
ERROR 



True 



Select the NSTAT=65 solution 
from Phase 1 with the highest 
\alue for the objecme 
function 



2 1 Run the NLP on this 
solution with all weights 
that are zero clamped to 
zero. 



NSTATot 



REPORT 
ERROR 



NSTAT =1 

© 



FIG. 7 







True 


i »■ 

I All weights > minimum 7 

I 

i . . 


► 



•0 



False 



3 0. Does each asset class hd\c at least 
one security meermg the minimum initial 
investment requirement 0 



False 



True 



3 1 Set all nonzero weights 
in current opumai solution to 
minimum and rerun 



j 3 2 Clamp to zero ail 
weights in current optimal 
solution that do not meet 
minimum investment. 



If both solutions fat!, repeat this 
branch with problem extended to 
allow asset class bounds to be 
exceeded (as in 3.4) 

Otherwise, select best feasible 

<;n!nf!on and on rn nha<sr> 4 



FIG. 8A 



FIG. 8B 
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At this point a solution has been identified where at least one asset class 
has zero securities meeting the minimum inv estment requirement Proceed 
as follows 

For asset classes vvith at least one security meeting its minimum 
requirement, include ail securities meeting the min, and limit their lower 
bound to the mm, For ail assets in these classes not meeting the mm, set 
the upper bound to zero 

For asset classes with no securities meeting the mm. identifv the secunt\ 
uitb the largest weight and set it^ lower bound to the nun investment 

Sohe the redefined NLP 



\ST\T 




For each asset class (excluding "other"), add two variables representing 
positive and negative deviations (both will take on positive values between 0 
and 1). Thus in the current formulation, add 12 variables to the basic 
problem. In addition to using these to force the lhs to meet the rhs of the 
asset class constraints, they are appended to the objective function with a 
penalty determined by the fixed parameter Gamma (currently set to 100). 
Thus, the sum of deviations is minimized. Note that negative deviations 
from the cash target are currently penalized at 2\. At some point, we should 
specify vvhat the maximum allow able deviations for each asset class 

Solve the redefined NLP 



NSTAT 



NSTAT <> i 



REPORT ERROR 



Phase 4 



At this point an optimal solution has been identified from the NLP. This cleanup phase 
vwfl make the transformations necessary to link the results with the finai user interface. 



REPORT ERROR I 



r 



There are three identifiable eases to be addressed \\ here a solution can not be reported to 
the u^er, 

1) The NLP can not soKe due to the lack of funds — it is possible, a prion, to determine 
when there 1^ not a feasible solution due to the fund mimmums and how the 
discretionary' funds are allocated across the tax^nontax accounts. A reporting 
mechanism should be designed to account for this case. An alternative to consider is 
that for these cases ue recommend an asset allocation fund appropriate for the user's 
risk tolerance, such as the Schwab MarketTrack funds. 

2) The NLP does not indicate a solution in any phase before the minimum investment 
criteria are considered — we do not anticipate that this will occur, given the 
robustness of the iterative solution method, but it cannot be ruled out. Thus an error 
reporting mechanism must provide for the case where the NLPjust does not work. 

3) There are situations in phase 3 where it is possible for the NLP to fail. Although the 
flexibility of the 3.4 technique should preclude most alt of these, an error is still 
possible. We need an error reporting mechanism for this possibility. In these cases, 
however, we may be able to modify the solution approach to resolve the error. 

The error reporting system will include a log file where the step at which the solution 
failed is identified This will help us identify errors where a change to the algorithm can 
be made to avoid similar errors in the future 
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